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ABSTRACT 


V\ 

'"'This  work  addresses  the  design  and  analysis  of  the 
Pitch,  Yaw  and  Roll  autopilot  for  application  to  the 
Bank-To-Turn  (BTT)  missiles.  At  first,  the  linear 
uncoupled  channels  were  designed  and  analyzed  according  to 
the  desired  requirements.  Utilizing  the  uncoupled 
channels,  the  linear  coupled  autopilots  were  designed,  not 
including  the  inertial,  kinematic  and  aerodynamic 
cross-coupling.  Then,  the  nonlinear  CBTT  autopilots  were 
designed  and  analyzed,  using  the  linear  CBTT  (coordinated 
Bank-To-Turn)  models,  which  now  have  coupled  with 
kinematic,  inertial  and  aerodynamic  cross-coupling.  The 
minimization  of  the  above  kinematic  and  inertial  coupling 
and  their  effects  were  completed  using  feedbacks  of 
angle-of-attack  and  rate  of  angle-of-attack  in  the  Pitch 
autopilot. 
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TABLES  OP  SYMBOLS  AND  ABBREVIATION 

Bank-to-Turn 

coordinated  Bank-to-Turn,  minimum  sideslip, 
positive  a,  4e  180  deg 

coordinated  Roll-During-Turn 


rolling  moment  coefficient 

C  slope  of  curve  of  rolling  moment  coefficient, 

•  ct  vs  • 

C  change  in  C,  per  degree  roll  control  incidence, 

*  •  •  * 


change  in  C  per  degree  yaw  control  incidence, 
*Y 

pitching  moment  coefficient 

slope  of  curve  of  pitching  moment  coefficient 
C  vs  a 

m 

change  in  Cm  per  degree  pitch  control 
incidence,  sp 

normal  force  coefficient 

slope  of  curve  of  normal  force  coefficient  CN 
vs  a 

change  in  CN  per  degree  pitch  control 
incidence,  6p 


CR  yawing  moment  coefficient 

Cn  slope  of  curve  of  yawing  moment  coefficient,  C_ 

e  vs  b 


pitch  angular  acceleration  about  YB 
constant  or  equilibrium  pitch  angular  rate 
yaw  angular  rate  about  Z0 

yaw  angular  rate  command  (coordination  command) 
yaw  angular  acceleration  about  ZB 

2 

reference  area  for  coefficients  *  »  ft 

Skid-to-Turn,  roll  attitude  stabilized 

velocity  component  in  xB  direction 

velocity  component  in  Y0  direction,  assumed  to 
be  constant 

constant  missile  flight  path  velocity 

missile  velocity  vector 

velocity  component  in  ZB  direction 

body-fixed  roll  axis,  along  axis  of  symmetry, 
positive  forward 

body-fixed  pitch  axis,  positive  starboard 

vehicle  axis  in  local  horizontal  direction, 
approximated  as  inertial  axis 

body-fixed  yaw  axis,  forms  right  handed 
orthogonal  system  with  XB  and  YB 

vehicle  axis  in  downward  direction  along  local 
gravity  vector,  approximated  as  inertial  axis 

achieved  normal  acceleration  in  z v  direction 
commanded  normal  acceleration  in  ZB  direction 


achieved  normal  acceleration  in  YB  direction 


achieved  normal  acceleration  in  zD  direction 

D 
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C  change  in  Cn  per  degree  yaw  control  incidence, 

«y  «Y 

C  change  in  C  per  degree  yaw  control  incidence, 

"4R  *R 

Cy  side  force  coefficient 

C„  slope  of  curve  of  side  force  coefficient  Cy  vs 

Ye  S 

Cy  change  in  Cy  per  degree  yaw  control  incidence, 

0  y  0  y 


change  in  Cy  per  degree  roll  control  incidence, 
aR 


d 

1YY 

Izz 

*xx 

KA 

Kyp 

LBTT 

P 

P 

Pe 

POC 

q 

q 


reference  length  for  cefficients  »  2  ft. 

moment  of  inertia  about  Y0  axis 

moment  of  inertia  about  ZB  **is 

moment  of  inertia  about  Xg  axis 

autopilot  pitch  acceleration  error  gain 

CBTT  autopilot  coordination  branch  gain 

limited  Bank-to-Turn  may  or  may  not  be 
coordinated,  positive  and  negative  a,  90 

or  45  degrees 

roll  rate  about 

roll  acceleration  about  XB 

constant  or  equilbrium  roll  angular  rate 

preferred  orientation  control 
dynamic  pressure 
pitch  rate  about  Yb 


achieved  normal  acceleration  in  Yv  direction 

normal  acceleration  command  from  guidance 
computer  in  Zy  direction  plus  anti-gravity  bias 
command 

roll  attitude  command  from  guidance  computer, 
zero  degrees  in  Zv  direction  and  90  degrees  in 
Yv  direction 

roll  attitude,  zero  degrees  in  -Zy  direction 
and  90  degrees  in  Yy  direction 

Elevation  Euler  Angle,  second  rotation, 

(q  cos#  -  r  sin#)  dt 

pitch  control  incidence  (positive  tail 
incidence  produces  negative  pitching  moment)* 

commanded  pitch  control  incidence,  sp 


yaw  control  incidence  (positive  tail  incidence 
produces  negative  yawing  moments)* 

commanded  yaw  control  incidence,  «Y 

roll  conntrol  incidence  (positive  tail 
incidence  produces  positive  rolling  moment) 

commanded  roll  control  incidence,  «R 


constant  or  equilibrium  angle-of-attack 

angle-of-attack 

ange-of-attack  rate 

modified  form  of  estimated  angle-of-attack  for 
autopilot  coordination  command 

angle  of  sideslip 

sideslip  angular  rate 
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I.  INTRODUCTION 


In  recent  years,  the  application  of  Bank-to-turn 
guidance  to  missiles  has  generated  considerable  interest. 
This  interest  has  been  motivated  by  the  drag  reduction, 
which  the  Bank-to-turn  (BTT)  controlled  missiles  offer, 
over  the  conventional  cruciform  roll  stabilized, 
skid-to-turn  controlled  missiles. 

Bank-to-turn  steering  may  provide  improved  performance 
for  missile  systems.  In  addition,  the  ramjet  engine  inlet 
configuration  is  believed  to  provide  a  greater  range 
capability  than  other  inlet  configurations.  BTT  control  is 
required  to  satisfy  the  sideslip  constraints  imposed  by 
chin  inlets.  Despite  these  facts,  there  are  unanswered 
questions  concerning  BTT  systems  stability  during  homing, 
guidance  performance,  autopilot  and  guidance  logic  design 
and  subsystem  requirements  which  must  be  investigated 
before  Bank-to-turn  steering  can  be  considered  a  viable 
method  of  control  for  high  performance  missiles. 

Many  missile  programs  were  initiated  during  the  past 
decade  to  improve  their  capability  via  Bank-to-turn 
control.  The  results  have  advanced  the  understanding  of 
the  different  missile  subsystems.  In  the  autopilot  area, 
many  types  of  autopilots  have  been  found  which  force  the 
missile  to  roll  or  bank  so  that  the  steering  maneuver 
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occurs  with  the  airframe  oriented  in  a  specified  or 
preferred  direction  with  respect  to  the  incoming  airstream. 
This  entire  class  of  autopilots  may  be  referred  to  as 
preferred  orientation  control  autopilots. 

Also,  simplified  guidance  studies  which  neglect  radome 
effects  and  assume  coordinated  missile  motion  have  shown 
that  CBTT  can  provide  acceptable  performance  with  roll 
rates  that  are  not  excessive  for  autopilot  control.  These 
studies  were  made  for  a  medium  range  area  defense  mission 
and  a  long  range  suppression  mission  and  considered  both 
high  lift  and  moderate  lift  configurations  [Ref.  1]. 

The  present  work  reviews  the  design  and  analysis  of 
the  pitch.  Yaw  and  Roll  autopilot  for  application  to  the 
Bank-to-turn  (BTT)  missiles  developed  in  [Ref.  2].  At 
first,  the  linear  model  of  the  uncoupled  autopilot  channels 
will  be  designed  and  analysed  in  the  critical  areas  o£ 
concern  regarding  the  response  of  CBTT  control  for  the 
circular  and  elliptical  airframes.  Then,  the  designed 
autopilot  channels  will  be  coupled  using  suitable 
aerodynamic  models  and  a  measure  of  sideslip  control  is 
obtained  by  roll  command  to  the  linearized  CBTT  autopilot. 
The  relative  stability  of  the  autopilot  branches  and  means 
for  improving  stability  will  be  discussed.  In  a  later 
chapter  the  non-linear  three  dimensional  model  will  be 
designed,  using  a  three  dimensional  aerodynamic  model  and 


the  same  control  laws  described  in  the  linear  studies, 
except  for  a  minor  modification  to  the  coordinating  branch 
dependence  on  Angle-of-Attack  and  also  the  inclusion  of 
Anti-Gravity  Bias.  Next,  the  kinematic  and  inertial 
coupling  effects  between  pitch  and  yaw  dynamics  will  be 
studied.  Transients  in  maneuver  plane  acceleration  are 
caused  in  the  above  cross-coupling  effects  and  are  in  the 
form  of  overshoots  and  undershoots.  Coupling  transients 
will  be  reduced  by  increasing  pitch  stability.  The 
technique  of  feedbacks  of  Angle-of-Attack  and  rate  of 
Angle-of-Attack  for  the  pitch  autopilot,  will  be 
considered,  in  order  to  decrease  kinematic  and  inertial 
cross-coupling  effects.  Finally  conclusions  and 
recommendations  for  future  study  will  be  stated. 


II 


LINEAR  DESIGN  AND  ANALYSIS  OF  UNCOUPLED  AUTOPILOTS 


FOR  CIRCULAR  AND  ELLIPTICAL AIRFRAME 


A.  GENERAL 

The  linear  design  and  analysis  started  with  uncoupled 
autopilot  channels,  i.e.  aerodynamic,  kinematic  and  control 
law  cross-coupling  between  Pitch,  Yaw  and  Roll  channels 
were  removed.  The  autopilot  design  technique  was 
classical,  using  a  combination  of  frequency  response  and 
root  locus  techniques  [Ref.  2],  providing  the  range  of 
required  missile  body  angular  rates,  control  motions  and 
minimizing  the  influence  of  aerodynamic  variations  on 
desired  responses.  To  achieve  the  desired  maneuver  plane 
acceleration  response  for  the  CBTT  autopilot,  both  Pitch 
and  Roll  uncoupled  autopilots  were  designed  to  have  a  time 
constant  of  0.5  sec  in  order  to  achieve  the  desired 
maneuver  plane  speed  of  response.  The  Yaw  channel  which 
follows  the  Roll  channel  motion  to  produce  desired 
coordination  (i.e.  minimization  of  sideslip  angles) ,  was 
designed  to  have  a  more  rapid  response  than  the  Roll 
channel,  having  0.39  sec  for  circular  and  0.36  sec  for 
elliptical  time  constant  in  the  maneuver  plane  acceleration 
response.  The  whole  design  was  done  for  the  constant 
flight  condition  of  60  KFT  altitude  and  mach  number  3.95. 
The  design  discussed  here  was  developed  by  Arrow  [Ref.  2]. 
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In  this  chapter  the  airframe  configurations,  the 
aerodynamic  models  and  control  laws  for  the  Pitch,  Yaw  and 
Roll  channel  are  reviewed.  The  transfer  functions  for  both 
circular  and  elliptical  configuration  are  derived.  An 
analysis  for  every  channel  in  regard  to  the  acceleration 
responses,  the  body  angular  rates  and  control  surface 
deflection  is  made.  Also,  a  comparison  of  airframe 
responses  for  every  channel  is  made,  considering  the 
relative  stability  and  required  control  surfaces. 

A  general  block  diagram,  of  the  channels  and  the 
command,  which  are  applied  to  the  control  laws  is  shown  in 
Figure  2.1.  A  number  of  the  Figures  shown  in  this  report 
have  been  taken  from  [Ref.  2]. 

B.  AIRFRAME  CONFIGURATIONS 

The  U.S.  National  Aeronautics  and  Space  Administration 
maintains  a  continuing  research  effort  in  missile  related 
aerodynamics.  In  this  effort  monoplanar  missile  concepts 
have  been  considered,  in  part,  because  of  the  relatively 
low  geometric  profile  that  can  be  achieved  for  convenient 
carriage  [Ref.  1].  The  two  Airframe  Configurations  studied 
in  this  research  were  taken  from  [Ref.  1]  and  are  shown  in 
Figure  2.2  and  Figure  2.3.  The  circular  cross  sectional 
body  has  a  closure  ratio  Abase/Amax  0.69  and  the  Amax 
occurs  at  68.0%  body.  The  Elliptical  has  an  exact  3:1 
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Autopilot  (BTT) 


Configuration  -  1/6  Scale 


Model  of  Elliptical  Configuration  -  1/6  Scale 


elliptical  cross  section,  and  its  cross  section  area  is  the 
same  as  the  area  of  circular  body.  The  span  of  both 
airframes  is  the  same.  The  wing  area  and  span  for  the 
circular  concept  were  chosen  as  typical  of  current 
maneuvering  missiles.  The  wing  for  the  elliptical 
configuration  was  determined  by  projecting  the  elliptical 
body  on  the  circular  bodywing  planform.  The  resultant 
exposed  wing  planform  then  became  the  wing  for  the 
elliptical  body.  The  tails  for  each  design  had  a  dihedral 
of  30°  degrees  as  shown  in  Figure  2.2.  In  order  for  tail 
deflection  to  be  compatible  with  complex  surfaces  of  the 
afterbody  of  the  elliptical  configuration,  the  tail  hinge 
line  was  skewed  such  that  a  10  degrees  deflection  measured 
at  the  body-tail  juncture  had  a  resultant  7.04  degrees 
surface  deflection,  as  shown  in  Figure  2.3  and  [Ref.  1]. 
Comparison  of  the  elliptical  cross  section  concept  with  the 
circular  cross  section  concept  indicates  the  following: 

1.  About  25  percent  more  normal  force,  that  is  nearly 
independent  of  angle  of  attack  can  be  achieved  at 
supersonic  speeds. 

2.  Values  of  the  longitudinal  stability  parameter  Cm 

“a 

are  more  positive,  with  more  pronounced  nonlinearities  in 
pitching  moment  occur ing  at  subsonic  speeds. 

3.  Levels  of  directional  stability  are  increased  and 
are  more  compatible  with  levels  of  longitudinal  stability. 


4.  Noticeably  more  yaw  control  is  available,  although 
suitable  locations  for  tails  on  the  body  are  more  limited 
because  of  the  geometry  of  the  elliptical  body. 

C.  UNCOUPLED  LINEAR  PITCH  CHANNEL 
1.  Aerodynamic  Model 

The  initial  phase  in  the  design  of  the  CBTT 
autopilots  involves  the  design  of  the  uncoupled  pitch 
channel.  The  linear  pitch  channel  dynamic  model  which  is 
used,  is  shown  in  Figure  2.5. 

This  figure  is  taken  from  [Ref.  2].  The  above 
model  is  linearized  for  stability  studies.  The  following 
three  assumptions  were  made: 

a.  Plane  X0-Zg  is  the  maneuver  plane,  shown  in 
figure  2.4 

b.  Missile  is  trimmed  in  pitch,  that  is  My=0,  at 
fixed  values  of  «,  q  and  *p. 

c.  Missile  roll  rate  is  constant.  Using  the 
foregoing  assumptions,  and  the  approximation  that  the 
variation  in  forward  velocity  (u) ,  is  equal  to  zero,  a 
simplified  set  of  longitudinal  equations  of  motion  are 
constructed  applicable  to  short  period  motion.  These 
equations  describe  the  two  degree  of  freedom,  short  period 
mode  longitudinal  motion. 

Linearized  aerodynamic  derivatives  for  M*3.95 


are  provided  in  Table  I. 
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Aerodynamic  Sign  Convection 


q(  rad/sec) 


TABLE  I 


Linearized  Aerodynamic  Derivatives  (M  *  3.95) 


Circular  Elliptical 


a  =  0° 

a“ 

10° 

a  * 

20° 

a  ■ 

>  0° 

a  =* 

10° 

a  * 

20° 

-  .065 

- 

.082 

- 

.111 

- 

.043 

- 

.054 

- 

.064 

s 

-  .025 

- 

.019 

- 

.003 

+ 

.024 

+ 

.024 

+ 

.032 

s 

0 

- 

.009 

- 

.020 

0 

- 

.027 

- 

.040 

CYS 

4V 

+  .021 

+ 

.022 

.028 

+ 

.016 

+ 

.015 

+ 

.019 

s 

SY 

-  .050 

— 

.053 

.062 

— 

.042 

— 

.039 

— 

.045 

s 

5V 

0 

— 

.016 

— 

.038 

0 

— 

.010 

.023 

CY. 

*R 

0 

— 

.009 

.022 

0 

— 

.006 

— 

.014 

Cn{ 

"R 

0 

+ 

.018 

+ 

.044 

0 

+ 

.014 

+ 

.032 

% 

®R 

+  .031 

+ 

.035 

+ 

.044 

+  . 

.023 

+ 

.023 

+ 

.029 

CN 

a 

+  .15 

+ 

.17 

+ 

.22 

+ 

.18 

+ 

.22 

+ 

.30 

\ 

+  .04 

+ 

.04 

+ 

.05 

+ 

.02 

+ 

.02 

+ 

.025 

Cm 

a 

-  .060 

— 

.065 

— 

.118 

+ 

.015 

+ 

.0137 

— 

.125 

Cm6 

®D 

-  .080 

- 

.095 

- 

.115 

- 

.055 

- 

.055 

- 

.075 

Reference  C.G.  at  0.6 
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.2.  Pitch  Control  Law 

The  pitch  control  laws  for  the  circular  and 

elliptical  airframes  are  shown  in  Figure  2.6  as  given  in 

the  [Ref.  2].  Lag -Leads  were  used  to  prevent  guidance 

noise  saturation  problems.  An  integrator  was  used  in  the 

acceleration  error  branch  of  the  circular  control  law  to 

satisfy  the  guidance  requirement  of  zero  steady  state 

error.  In  the  elliptical  control  law,  an  integrator  was 

placed  in  the  rate  error  path  (i.e  path  after  the  adder 

(q+y)  •  as  shown  in  Figure  2.6)  for  reduction  of  the 

acceleration  response  overshoot  below  10  percent  and  a  gain 

K*1 . 1053  was  placed  in  series  with  the  acceleration  command 

for  zero  steady  state  acceleration  error. 

A  normal  acceleration  command  nz  equal  to  1  gee, 

c 

is  applied  to  the  pitch  control  law  which  uses  measurements 

of  missile  body  pitch  angular  rate  (q)  and  pitch  normal 

acceleration  (nz)  to  determine  the  required  actuator 

command  (6_  ).  The  actuator  is  modeled  as  a  first  order  lag 
*c 

at  188.4  rad/ sec  and  is  shown  in  Figure  2.5. 

3.  Transfer  Functions  of  Aerodynamic  Models  for 
Circular  Airframe 

The  equations  which  are  represented  by  the  block 


diagram  of  the  aerodynamic  model,  of  the  uncoupled  Pitch 
channel  (Figure  2.5)  are: 


Itch  Control  Laws 


tv:* 


W  [$««] 

0L=  5T3[Cj  +  3a.Xt\^/ v]  [4e$] 


[foJ/5«cJ 
(II. C. 3-1) 

(II. C. 3-2) 
(II. C. 3-3) 


Utilizing  Laplace  Transformation,  the  above 


equation  becomes: 

$r  AM  (Ctng-a  -hS^C^-dp) 

IYy 

“9^  (51.3  Cw$p  +  C^-  a)/w 

5a=  5T.3(q  +  3iatl*/V) 


(II. C. 3-4) 

(II. C. 3-2) 
(II. C. 3-5) 


Introducing  equation  II. C. 3-2  to  II. C. 3-5  the  last 
equation  yields. 


3  $<L  (  +  5?*3  Cn»JP  *^P  ) 

tyy 


(II. C. 3-6) 


(II. C. 3-7) 

Rearranging  above  equations,  in  terms  of  q  and  a, 
we  obtain: 

S<\- (-^ -C^) • • «  =  ( ” C»S9 ) ’ S?  (II.C.3-8) 

=  (-  ^  ^  }*  UI.C.3-8, 

Applying  the  Cramer's  rule  in  the  system  of 

equation  II.C.3-8  and  II. C. 3-9  the  transfer  functions  q/5p 

and  a/ft  are  obtained. 

P 
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'■/■.•'•.•'■.•VvV 
^  »  -’■I .'-1  ^ 


.-k..--  -  v  v 
\  .v  - 


As  referred  to  in  Appendix  A  and  Table  I.  After 


substituting  the  aforementioned  values  into  equation 
II. C. 3-10,  it  becomes: 

-3 _ — z  S  -  O  -ISQ  f  4 ir  \ 

^P  22.5MS  XHd*Sa+ 1.1615  XiO^S  +i  '  (II. C. 3-12) 

Introducing  equation  II. C. 3-12  into  equation 


II. C. 3-2  the  last  equation  yields: 


^ - 7s  ~ 


3283.23  (QSfij-CmS,- 

Iyy  •  W 


t* 


ST-3  Cf  jfd  Gna 
Iyy 


5T-5q^ 

W 


(II. C. 3-13) 

Substituting  the  above  mentioned  values  into 
equation  II. C. 3-13  and  minor  manipulating,  it  yields: 


ft  -  iOS.OmjdV  -  19.$2 

*P  ~  22.  SHSXiCT5-  S*+3.3 C35X10’5  +  i 


(%<ts/ rad) 

(II. C. 3-14) 


The  equations  II. C. 3-12  and  II. C. 3-14  are  the 
pitch  aerodynamic  transfer  functions  for  the  pitch  angular 
rate  q  about  the  and  the  achieved  maneuver  acceleration 
in  the  zB  direction  (nz) . 

4.  Design  Approach  and  Analysis  For  Circular  Airframe 
A  Fixed  Flight  condition  was  selected  for  these 
preliminary  performance  studies  of  circular  and  elliptical 
airframes.  The  selected  flight  condition,  at  60KFT 
altitude  and  mach  3.95,  provides  a  sufficiently  low  dynamic 


pressure,  so  that  missile  maneuvers  will  result  in  large 
enough  angles-of-attack,  to  exercise  side  slip  control. 

The  above  fixed  condition  was  selected,  in  order  to  reduce 
the  complexities  of  design  of  the  CBTT  autopilots. 

The  uncoupled  autopilot  design  technique  was 
classical,  using  a  combination  of  frequency  response  and 
root  locus  technique  [Ref.  2]  to  achieve  practical 
bandwidths  (i.e.,  sufficient  high  frequency  attenuation) 
and  therefore  provides  the  range  of  required  missile  body 
angular  rates  and  control  motions.  The  analysis  of  the 
CBTT  autopilot  channels  is  based,  on  the  time  responses  of 
manuever  plane  acceleration,  body  angular  rates  and  tail 
incidence  angles. 

For  purposes  of  analysis,  a  CSMP  program  was 
written  (Appendix  D)  using  the  equations  which  are 
represented  by  the  block  diagrams  of  aerodynamic  model  and 
linear  uncoupled  pitch  control  law,  as  are  shown  in  figures 
2.4  and  2.6. 

The  equations  of  the  uncoupled  linear  pitch 
aerodynamic  model  and  control  law  are: 
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x  =  -150  X+  450^ 


(II. C. 4-2) 


(2)  Acceleration  Compensator  Equation 

-««•  fe*)  ,A 

T  —  - - — - - r-'l^-X)  (II. c.  4-3) 

r  f  .  »  A 


H- 


Utilizing  inverse  Laplace  Transformation  and 


rearranging  the  above  equation,  it  becomes: 

V+5V  =  -1.3534  C^2c-X)  -mSXJO^ic-x) 


(II. C. 4-4) 


Using  state-space  representation  of  a  system, 
in  which  the  forcing  function  involves  derivative  terms 


[Ref.  5:  pp.  675-678]  one  obtains: 


Y|  -13531  (K^-x) 

\  =-5Y2  +6.532  •  (Hfc-X) 

Y  =  \ 

(3)  Rate  Compensator  Equation 


(II. C. 4-5) 
(II. C. 4-6) 
(II. C. 4-7) 


^Pc  = 


(~ohi+i) 


(V*<) 


(II. C. 4-8) 


yields : 


Using  Inverse  Laplace  Transformation,  it 

ip,  r  -<M43-Spc-  3M0.5SX40  )  -  2.230S  •  (Y4-X<) 

rc  *•  (II.C.4-8a) 


Actuator  Equation  (Figure  2.5): 

■= - — — -  -Jpc 

J  I  I 


(II. C. 4-9) 


+  1 
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(II. C. 4-10) 


<5>p  =:  -iW.4*  <5p  +  188.4  Sfc 

c.  Equations  of  Aerodynamic  Model 

In  order  to  achieve  the  plots  of  the  time 
response  of  q,  nz,  6  in  the  linear  uncoupled  pitch 
channel ,  the  transfer  functions  which  are  represented  by 
the  aerodynamic  model  and  derived  in  section  II. C. 3,  were 
used : 

9  -1.33  CIS  -0.169  (II.  C.  4-11) 


25.5i|5Xl6,52  +  3.3635X-l63S  + 1 


or 


<1  (  25.545  X 10-5*+  3.3633  X 1 6  •  S +1 )  =  (-  1.33CJ -5  -  0.159>  Sp 

(II.C. 4-lla) 

Rearranging  and  minor  manipulating  equation 
II.C. 4-lla  it  yields: 

<$  +  1Mg.imi03q  +  M<(.355C  9  =  59.2635  £p  -*.052  £p  (n.c.4-12) 

Utilizing  state  representation  of  a  system  in 

which  the  forcing  function  involves  derivative  terms,  as  in 

section  II.C. 4,  one  obtains: 


Xi  -  X2  -  59.2*35-^ 


(II.C. 4-13) 


-3 


jL  -  “44.3556  X^- 449. 182X10  Xa+J.m£p  (II.C. 4-14) 


(II.C. 4-15) 


q  = 

Also,  rearranging  and  manipulating  in  the  same 
way,  equation  II.C. 3-12  yields: 

1*  -  12  -  J0J.Q01  XIO5.  5p  (II.C. 4-16) 


^2  = -^.3354  -i.WSSXIO^p  (II.C. 4-17) 


-3 


»  ■  ■  <• 
\>v 


^2  -  "  (*1  +  M.?05  <Tp  ^  (II. C.  4-18) 

Utilizing  state-space  representation,  all  the 


aforementioned  equations  can  be  modeled  in  the  following 
ninth  order  system: 
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0  . 

• 

Vi 

• 
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0  0 
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-5 
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Y2 
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0344 

0 

0  0 

-0.461 
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* 

L 
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(II. C. 4-19) 


where  the 


The  state 
figure  2.6 


state  variables  are: 

,  X2:  pitch  angular  rate  state  variables 
(section  II. C. 4) 

Y^,  Tj!  Achieved  Body  Acceleration  state 
variables  (section  II. C. 4) 

:  output  of  Acceleration  filter 


X 

6. 


input  command  in  actuator  network 
tail  incidence  angle 


X,  Ylf  5p  of  control  law  are  shown  in 
c 
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Executing  the  program  of  Appendix  D,  using  an 
input  step  function  which  represents  the  "lGee  Command"  and 
trim  angle  of  attack  ae  =0,  the  time  response  plots  of 
pitch  normal  acceleration  nz,  pitch  angular  rate  q  and 
pitch  tail  incidence  «p  for  circular  airframe  are  obtained. 

Figure  2.7  shows  the  pitch  acceleration  response 
of  the  circular  airframe  due  to  a  one  Gee  acceleration 
command,  when  the  Missile  Aerodynamic  Models  are  linearized 
about  zero  angle  of  attack  (a).  The  response  has  a  0.5 
seconds  time  constant,  steady  state  error  0.005  and  2.5% 
overshoot,  which  are  according  to  the  requirements  referred 
to  in  appendix  C. 

Figures  2.8  and  2.9  show  the  required  body 
angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response.  These  responses  matched  those 


presented  in  [Ref  2]. 

5.  Design  Approach  and  Analysis  For  Elliptical 
Airframe 


The  design  technique  used  for  the  elliptical 


airframe,  is  the  same  as  that  of  the  circular  airframe. 


The  analysis  of  the  Pitch  channel  for  the  elliptical 


airframe  was  based  again,  on  the  time  response  plots. 


For  the  above  purposes,  a  CSMP  program  was 
written,  as  in  the  case  of  circular  airframe,  referred  to 
in  Appendix  D. 


44 


MtJi  Ul  dJI  1«  IJI  IJI  }«  U  I4D  3.n  3.00  3J* 

TCSBCJ 

Fig.  2.7  Pitch  Normal  Acceleration  (nz)  vs  Time  (t) 
Uncoupled  Pitch  Channel 
Circular  airframe  (lGee  Command,  a  =0). 


MMi  M  Ml  IJ>  IJi  IJI  IJI  M  U  3.M  3.7S  3.M  3.S 

Ttsec) 

Pig  2.8  Pitch  Angular  Rate  (q)  vs  Time  (t) 
Uncoupled  Pitch  Channel 
Circular  airframe  (lGee  Command,  a  =0) 


The  equations  of  the  uncoupled  linear  Pitch 
channel  for  elliptical  airframe ,  which  are  used  in  the 
above  CSMP  program,  are: 

a.  Equation  of  Control  Law 

The  control  law  of  the  uncoupled  linear  Pitch 
channel  for  the  elliptical  airframe  is  very  similar  to  the 
control  law  for  the  circular  airframe,  with  only 
differences  in  the  gain  of  acceleration  and  rate 
compensators.  Elliptical  control  law  is  shown  in  figure 
2.6.,  taken  from  [Ref.  2]. 

(1)  Acceleration  Filter  Equation.  The 
Acceleration  filter  equation  is  the  same  as  in  section 
II. C. 4 


(2)  Acceleration  Compensator  Equation 

Y=  (-O-OOS/^+lV  ( 1/1053 vx)  (II. c. 5-1) 

Utilizing,  inverse  Laplace  Transformation 
and  rearranging  II.C.5-1  it  becomes: 

530.5^  XIoNlic  +0.MS.X  (ii.c.5-2) 

(3)  Compensator  Equation 

5ft  -  -5.W-  (■»■**) 

V-X,  '  £ 

(II. C. 5-3) 

Using  inverse  Laplace  Transformation,  it 


yields: 

Spc=  -3M.?5X*03(y-X|)  -3.0?(Y-X<) 


(II.C.5-4) 


b.  Actuator  Equation  (Figure  2.5) . 


As  in  section  II.C.4,b 

c.  Equation  of  Uncoupled  Linear  Pitch 
Aerodynamic  Model 

The  equation  for  pitch  angular  rate  was 
derived,  as  in  section  II. C. 3,  the  only  differences  are  in 
the  constants  W=2475(lbs), 

Iyy=790  (Slug-ft2)  ,  C  =0.015,  C  =-0.055,  C  =0.02  and 

Ul  Ut  .  O  x 


C  =0.18,  as  referred  in  the  Appendix  A  and  Table  I. 
a 

Manipulating  according  to  the  inverse  Laplace 
Transformation  rules  and  rearranging,  one  abtains. 

~~  =  ^6frS+0.6tt)/(-W.6i4X4oV+  4595  X-fd3S  +0  (II. c.  5-5) 


Using  state-space  representation  of  a  system, 
in  which  forcing  function  involves  derivative  terms  [Ref. 


5:  pp.  675-678.]  It  yields: 

X*  =  X* -4i.29-Sp  (ii. c. 5-7) 

X2  =  ii.Mis-xl+4fg.o<ix^65x2-iMg^<J'p  (ii.c.5-8) 

S  =  Xi  (II.C.5-8a) 

Also,  Manipulating  with  same  way  the  equation: 

Hi  -  -211.4  H  XI 6*  5*+  31.6 _  (II. c.  5-9) 

Jp  ’  -8S.WXI03 sV  15*95X40*  S  +  i 


Yielding 

2t  S  l2  +528.3X10.  (II.  C.  5-10) 

h  =  +  (ll.c.5-11) 

y\ J  r  -  +  2-936  Sf)  (II. c. 5-12) 


.*• 

-■V 

-.  V 


1 

V 


*5* 


>£? 


Utilizing  state-space  representation,  all  of 
the  above  equation  can  be  modeled  in  the  following 
eighth-order  system: 


11.18  0.18 

O  O 

O  O 


O  O  o  1  O  0  0  0.5282 

0  O  11.28  0.18  0  0  0  -W.3 

0  0  150  0  H50  0  0 

0  0  0  0  0.M8  -4  0  0 


5 

Xi 

0 

*2. 

0 

2i 

0 

• 

iz 

X 

0 

y 

-0.53 

fa 

-aaos 

1 

Sp 

0 

(II. C. 5-13) 

where  the  state  variables  are: 

Xj_,  X2:  pitch  angular  rate  state  variables 


(section  II. C. 4) 


Y,,  Y~:  Achieved  Body  Acceleration  state 


variables  (section  II. C. 4) 
output  of  Acceleration  filter 


input  command  in  actuator  network 
tail  incidence  angle 


The  state  X,  X^,  Y^,  of  control  law  are  shown  in 

c 


Figure  2.6 
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Utilizing,  the  above  state  variables  system,  a 
CSMP  program  was  written  (Appendix  D)  to  achieve  the  time 
response  plots  of  nz,  q  and  6^.  As  input  a  step  function 
is  used,  which  represents  a  one  gee  command.  The  trim 
angle  of  attack  ( ag)  is  equal  to  zero. 

Figure  2.10  shows  the  pitch  acceleration 
response  of  the  elliptical  airframe  due  to  the  command.  The 
response  has  a  time  constant  of  about  0.5,  steady  state 
equal  to  zero,  and  overshoot  o%,  keeping  the  requirements 
referred  to  in  Appendix  C. 

Figures  2.  Hand  2.12  show  the  required  body 
angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response. 

• Comparison  of  Airframe  Response 

Figures  2.7  and  2.10  show  that  the  Pitch 
acceleration  response  of  the  circular  and  elliptical 
airframe  are  approximately  the  same  having  a  0.5  second 
time  constant  and  negligible  overshoots. 

Figures  2.8,2.9,2.11  and  2.12  show  that  to  achieve 
the  acceleration  response  the  elliptical  airframe  requires 
less  body  angular  rate  and  control  surface  deflection.  The 
elliptical  airframe  requires  less  control  surface 
incidences,  because  the  transfer  functions  have  a  higher  DC 
gain  and  poles  which  are  located  at  a  lower  frequency  due 
mainly  to  a  more  neutrally  stable  airframe. 
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Fig.  2.10  Pitch  Normal  Acceleration  (nz)  vs  Time  (t) 
Uncoupled  Pitch  Channel 
Elliptical  air£rame  (lGee  Command,  a  =0) 


TlSECI 


Fig.  2.11  Pitch  anc-iar  rate  (q)  vs  Time  (t) 
Uncoupled  titch  Channel 
Elliptical  airframe  (lGee  Command, 


DY(RRO) 


T(SEC) 


Fig.  2.12  Pitch  tail  incidence  (6  )  vs  Time  (t) 
Uncoupled  Pitch  Channel'3 
Elliptical  airframe  (lGee  Command,  a  =0) 


)]  is  directly 


The  DC  gain  HqS/w).(CN  .Cm  /  Cm  -  CN 

a  6  a  5 

P  P 

proportional  to  pitch  control  moment  C  and  inversely 

ra  5 

P 

proportional  to  magnitude  of  stability  margin  in  Pitch 

(i.e.,  C  /C-.  or  distance  from  center  of  pressure  to 
m  n 
a  a 

center  of  gravity) .  Thus  the  nearly  neutrally  stable 

elliptical  airframe  is  expected  to  have  a  higher  gain  than 

the  stable  circular  airframe  even  though  its  control  moment 

coefficient  is  slightly  smaller.  Higher  DC  gain  will 

require  less  control  surface  incidence.  The  zeros  of  nz/«p 

are  directly  proportional  to  C^  and  the  location  of  the 

a 

ratio  Cm  /  Cjj  or  the  distance  from  the  point  of 
6  6 
P  P 

action  of  tail  forces  to  the  center  of  gravity. 

Hence,  for  the  elliptical  airframe  which  has 

larger  C^  the  zeros  of  nz/5  are  located  at  a  higher 

frequency.  Lower  control  surface  incidences  and  higher 

frequency  n_/<5  zeros  of  the  elliptical  airframe  will 
z  p 

simplify  the  design  of  a  rapidly  responding  Pitch 
autopilot. 


D.  UNCOUPLED  LINEAR  YAW  CHANNEL 

The  purpose  of  the  Yaw  channel  autopilot  of  a  CBTT 
autopilot  is  to  minimize  sideslip  angle  or  provide 
coordinated  motion  between  Roll  and  Yaw  channels.  This  was 
accomplished  in  two  ways.  First,  the  uncoupled  Yaw  channel 
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autopilot,  (i.e,  Roll  and  Pitch  dynamic  effect  neglected) 
was  designed  as  a  regulator  (i.e  no  guidance  command  and 
with  rate  and  acceleration  feedback)  to  help  minimize 
sideslip  angle.  Second,  to  aid  in  sideslip  control,  the 
regulator  was  commanded  from  the  Roll  channel  as  explained 
in  chapter  III. 

1.  Aerodynamic  Model  Control  Law 

A  block  diagram  of  the  uncoupled  Yaw  channel  is 
shown  in  Figure  2.13,  taken  from  [Ref.  2],  In  this  diagram 
the  aerodynamic  model  and  the  Yaw  control  law  are  involved. 
Also,  the  Yaw  control  law  without  the  accelerator  is  shown 
in  Figure  2.14.  These  figures  were  abstracted  from 
[Ref.  2]. 

The  normal  acceleration  n  is  not  used  to 

yc 

command  the  CBTT  autopilot.  Instead,  it  is  used  for  the 
design  and  analysis  of  the  uncoupled  channel.  The  command 
used  by  the  coupled  system  is  shown  in  dashed  lines  and  is 
a  yaw  angular  rate  command,  r. 

w 

The  yaw  control  law  is  governed  by  missile  body 
yaw  angular  rate  (r)  and  yaw  normal  acceleration  n^.  The 
yaw  control  law  determines  the  required  command  6 

yc 

to  an  actuator,  which  is  approximated  as  a  first  order  lag 
at  30  Hz  (188.4  rad/sec). 

The  aerodynamics,  linearized  about  a  trim  ae 
angle-of-attack  equal  to  zero. 


Uncoupled  Yaw  Channel 


2.  Transfer  Functions  of  Aerodynamic  Model 

The  aerodynamic  characteristic  of  the  Yaw  channel 
are  modeled  by  a  second-order  system,  by  the  following 
equation,  where  the  states  are  r,  yaw  angle  and  6  , 


sideslip  angle: 

Tr  (q idlin')  •  (  (^-6 +  57*3  CtiSy Crcuilsec)  (xi.d.2-1) 

(dtj)  (ii. d. 2-2) 

and  also,  by  the  equation: 

Vly-  (q£/W)*  +5?.3Cy<yy.Sy  )  (3««S>  (II.  D.  2-3) 

substituting  equations  II. D. 2-3  into  II. D. 2-2,  the 
last  equation  becomes: 

8  z  ra  [-1+ y  (^f  (<V,-«  +  sn  c,viv)')  J  (ii .d. 2-4 ) 

Taking  the  Laplace  Transform  of  the  equations,  one 


(II. D. 2-3) 


obtains 

rs=-S 


•  ((*6  ■  8  4  5?.  3  C**  .<Jy  ) 


(II. D. 2-5) 


fc=  S*3  (-f  +  ^?  (2JL  fc,, -0+5*3  Cy;, .Si))) 


(II. D. 2-6) 


Rearranging  above  equations,  in  terms  of  q  and  a. 


they  become: 


rs  - 


d 


,C 


(II. D. 2-7) 


513 -v  +  (s-  ^  — %-h 


(II. D. 2-8) 
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Applying  the  Cramer's  rule  in  the  system  of  the 


equations  II. D. 2-7  and  II. D. 2-8  the  transfer  functions  r/6, 


and  B/6„  are  obtained. 

/  _ _ .•mJx’i  a 


r  _  y  ”**?*''■  C-^-C^C rir) 

17"  s7-  (jas. C,4).s  + 

1  ^  '  lit  /  (II. D. 2-9) 


*  "  *  (Jdl&L.).*- 


(II. D. 2-10) 


(57.3) . (32.3) 

K  «  -  *  0.48 

V 

q  *  1650  (lbs/ft2) 

S  »  3.14159  (ft2 


d  -2  (ft) 

Izz  *  810 (slug-ft2) 

w  -  2525  (lbs) 

C  »  -0.025 

B 
C 

n«y  -  -0.050 

Cy 

yB  -  -0.065 

Cy 

y«y  -  0.021 

a  «  0° 

e 

As  referred  to  in  Appendix  A  and  Table  I. 
Substituting  the  aforementioned  values  into  equation 
(II. D. 2-9),  it  becomes: 
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-5M-5MX-f6*52- 3.8 X 1 0'* S  4- 1  '  **c*  (II*D*2-11) 

Introducing  equation  II. D. 2-10  into  II. D. 2-3  and 

rearranging  the  last  one,  the  transfer  function  of  yaw 

normal  acceleration  in  terms  of  the  yaw  control  incidence, 

is  obtained: 


f  5?.SK(q*)*  ^  ,  , 

^ .  v  g[— Vy — IyyW 
*Y '  s,+pi£-.(((«ys_ 


s*sq^ 


(II. D. 2-12) 


After  substituting  the  values  into  the  equation 


(II. D. 2-12),  one  obtains: 

-13.W  S*+  1T.TC 
Jy  ”  + 1 


(%eti/tad) 


(II. D. 2-13) 


The  above  transfer  functions  are  used  for  finding 
the  time  responses  of  the  yaw  angular  rate  r  and  the  yaw 
normal  acceleration. 

Using  inverse  Laplace  transformation,  the  equation 


II. D. 2-11  and  II. D. 2-13  become: 


Y  + 51338X10* f  -  J83335T  s-36.76l  5y-j.i5  Sv  (il.D.2-14) 

My  +  51338  XlO**y  -183335  My  =  2.38*  Jy  -32S.G29$V  (II.  D.  2-15) 


The  equations  II. D. 2-14  and  II. D. 2-15  involve 
derivative  terms  in  the  forcing  function.  Using  rules  of 
state-space  representation  of  a  system  [Ref.  5:  pp 
675-678],  they  turn: 


61 


X<  =  X2-36.*6J-5v 

^2  r  18.335- X.j  -  5 1.338X10^2  "863.014  X40*iy 
V  -  Xi 

2,  z  ?2-  1225^  Xl03-<5v 

ii  =  1S.335-21  -51.333xl03*2-  2*J.$57-Sy 

Hv=  2<  +  2.3yf.£y 

3.  Equations  of  Yaw  Control  Law 

As  mentioned  before,  the  yaw  control  law  is 

governed  by  missile  body  yaw  angular  rate  (r)  and  yaw 

normal  acceleration  and  it  is  used  to  determine  the 

required  command  5y  to  the  actuator.  The  yaw  control  law 

c 

at  the  flight  condition  of  interest  (i.e,  Mach  3.95,  60kft) 

is  shown  in  Figure  2.14,  taken  from  [Ref.  2],  The  rate 

compensation  determines  the  high  frequency  attenuation  and 

is  used  to  minimize  aerodynamic  variations  on  the  quality 

of  regulation.  The  acceleration  compensation  determines 

the  acceleration  bandwidth  via  the  time  constant  of  the 

acceleration  response  due  to  a  step  command  of  acceleration 

at  nY  . 
c 

a.  Acceleration  Compensator  Equation: 

Y_ - &3m - (V+Hy)  (ii.d.3-1) 

0.2S  +  1  1  c  1J 

Inverse  transforming  and  rearranging  the 
II. D. 3-1,  it  turns  into: 

V  :  -5Y  +  i.59^3  (fly-'rtYc') 


(II. D. 2-16) 
(II. D. 2-17) 
(II. D. 2-18) 
(II. D. 2-19) 
(II. D. 2-20) 
(II. D. 2-21) 


(II. D. 3-2) 


VXV.s 


& 


Rate  Compensator  Equation: 

Svc=:  M‘*S  (y+r)  (II. D. 3-3) 


Utilizing  inverse  Laplace  Transformation 
and  introducing  equation  II. D. 2-16,  the  equation  II. D. 3-3 
turns  into : 

-  0.M85  (Y  +  Xi)  85  (Y+X-l)  (II.D.3-3a) 

c.  Actuator  Equation: 

<JV=  - i _  Jyc  (II. D.  3-4) 

- - -  +  1 

1S«.4 

Inverse  transforming  and  rearranging  the 


equation  II. D. 3-4,  it  turns  into: 

$Y  r-  188.4  •  5y  +  188*4  Svc  (II. D.  3-5) 

Using  state-space  representation,  all  the 
aerodynamic  model  and  control  laws  can  be  modeled  in  the 
following  seventh  order  system: 
x,  O  i  O  O  o  O  Xi  O 

X2  11-335  -<i5M  O  0  0  O  -0.W3  X*  O 

i,  0  O  O  1  0  0  -Q1226  2l  0 

il  -  0  0  11-335  -ft5iq  O  0  -J81.lt  •  2?  +  0  5|yc 

Y  O  0  159?  0  -5  0  3.8ft?  Y  -1.55? 


*2 

h 

k 

— 

Y 

ft 

(II. D. 3-6) 

where  the  state  variables  are: 

Xx,  X2:  pitch  angular  rate  state  variables 
(section  II. D. 2) 
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1 1  s 


Y1#  Y2:  yaw  acceleration  state  variables 
(section  II. D. 2) 

Y  :  output  of  acceleration  compensator 

6P 

c  :  input  command  in  the  actuator 

5p  :  yaw  tail  incidence 

The  state  X^,  Y,  6^  of  control  law  are  shown  in 

c 

Figure  2.14. 


4.  Design  Approach  and  Analysis  of  Circular  Airframe 

The  design  technique  used  for  the  Yaw  channel  is 
the  same  as  the  technique  used  in  the  Pitch  channel.  The 
analysis  of  it,  for  circular  airframe  is  based,  on  the  time 
responses  of  the  yaw  angular  rate  r,  the  yaw  normal 
acceleration  and  yaw  tail  incidence  6y. 

For  purposes  of  analysis,  a  CSMP  program  was 
written  (Appendix  E)  using  the  state-variable  system  which 
is  represented  by  the  block  diagram  of  the  yaw  aerodynamic 
model  and  control  law. 

Executing  the  referred  program,  using  as  input  a 
step  function,  which  represents  the  "one  Gee  Command"  and 
trim  angle  of  attack  ae»0,  the  time  response  plots  are 
obtained . 

Figure  2.15  shows  the  yaw  acceleration  response  of 
the  circular  airframe.  The  response  has  a  0.39  seconds 
time  constant,  7%  overshoot  and  steady  state  error  0.018 
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Fig.  2.15  Yaw  normal  acceleration  (n  )  vs  time  (t) 
Uncoupled  Yaw  Channel  y 

Elliptical  airframe  (lGee  Command,  a  =0) 
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which  are  according  to  the  requirements  referred  in 
Appendix  C. 


Figures  2.16  and  2.17  show  the  required  body  yaw 
angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response.  Again  these  results  matched  those 
of  [Ref.  2]  . 

5.  Aerodynamic  Transfer  Function  for  Elliptical 

Airframe 

The  transfer  functions  for  yaw  angular  rate  and 
yaw  normal  acceleration  for  the  elliptical  airframe,  were 
derived  in  the  same  way  as  in  section  II. D. 2.  The 
only  differences  are  in  the  following  constants 

w  =  2475  (lbs) 

Izz  »  853  (slug-ft2 

Cn  =  0.024 

n  B 

C  =  -0.042 

6 

Y 

CY  *  -0.043 
B 

cv  »  0.016 

6Y 

As  referred  to  in  Appendix  A  and  table  I. 
a.  Transfer  function  of  yaw  angular  rate  r. 

■j—  -  cn. D. 5-D 

O*  59.838X10  2.593X10 -S  +  i 


040  040  040  0.70  1.00  140  140  140  0.00  040  0.B  0.70  3.00  3.S 

TlSECJ 

Fig.  2.16  Yaw  angular  rate  (r)  vs  Time  (t) 
Uncoupled  Yaw  Channel 
Circular  airframe  (lGee  Command,  ae=0) 
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b.  Transfer  function  of  yaw  normal  acceleration  nY 

My  _  <UI.63SX<0-5-  \ 

£y  '  59.838X10'3- SJ4.2.593X3o’s+l  V  rad  )  (II-D-5"21 

Using  inverse  Laplace  transformation,  the  equation 
II. D. 5-1  and  II. D. 5-2  becomes: 

V+l»3333X40?f4*6.mM  =  -29.213.<5y-999.36X40-<fy  (II. D. 5-3) 
Yly  4M3.333X<03Tfiy  4&?iiMly  =  M8.6S3*Jy  (II. D.  5-4) 

Utilizing  state  representation  of  a  system  in 
which  the  forcing  functions  involve  derivative  terms,  as  in 
section  II. D. 2,  the  above  equations  yield: 

Xl  r  *2-  29.221  *<Jy  (II. D. 5-5) 

5(2=  -ium  X<  -^3.333Xto*X2 +26UW403.fy  (ii.d.5-6) 

1'-  Xl  (II. D.  5-7) 

i  =  £2-93.013X4  0*<5v  (ii.d.5-8) 

-43.333X1 6*  h -150-663? -Sy  (ii.d.5-9) 

-  £*4- i.9S1Voy  (ii. d. 5-10) 

The  aerodynamic  model  for  elliptical  airframe  is 
shown  in  Figure  2.13. 


6.  Yaw  Control  Equation  for  Elliptical  Airframe 


The  control  law  of  uncoupled  linear  Yaw  channel  for 


elliptical  airframe,  is  similar  to  the  control  law  for 
circular  airframe.  There  are  differences  in  the  gains  of 
the  acceleration  and  rate  compensators,  as  is  shown  in 


a.  Acceleration  Compensator  Equation 

y=  (  0.«3<35/o.25-S  +1)  •  (-H*  +  Tlv')  ( 1 1 .  D.  6-1 ) 

Utilizing  inverse  laplace  transformation  and 
rearranging  the  equation  II. D. 6-1,  it  turns  into: 

y  =  -4  V  +  135?M  •  (-nvc  +1M  (ii.d.6-2) 

b.  Rate  Compensator  Equation 

Svt:  *•<>«•  +  ( v-rl 

S  (II. D. 6-3) 

Inverse  Laplace  transforming  and  introducing 
the  equation  (II. D. 5-7)  into  the  equation  (II. D. 6-3),  it 
turns  into: 

<$YC  =  a608(y+X<)  +6.0$(Y+XA  (II. d. 6-4) 

c.  Actuator  Equation 

Same  as  the  equation  II. D. 3-5. 


Using  state-space  representation  all  equations  of  the 
aerodynamic  model  and  control  law,  can  be  modeled  in  the 
follow  seventh-order  system: 


x< 

0 

1 

O 

0 

0 

0 

-29.22 

*1 

0 

Xz 

-iU2 

-0.433 

0 

0 

0 

0 

0.26? 

X2 

0 

i 

O 

O 

0 

1 

0 

0 

-0*3 

h 

0 

% 

h 

r 

0 

O 

-0.433 

0 

0 

-150.6 

• 

h 

+ 

0 

* 

0 

0 

3.36 

0 

-4 

0 

6.553 

Y 

-3.36 

Stc 

609 

0.608 

2.04* 

0 

-1432 

0 

3.984 

Sk 

-104 

$y 

0 

0 

O 

0 

O 

1*8.4 

-1»4 

(II.D.6-55 
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where  the  state  variables  are: 


X1 '  x2: 


The  state  Xj.,  Y, 
2.14. 


yaw  angular  rate  state  variables 
(section  II. D. 5) 
yaw  normal  Acceleration  state 
variables  (section  II. D. 5) 

output  of  acceleration  compensator  network 
input  command  in  the  actuator 
yaw  tail  incidence 

Sy  of  control  law  are  shown  in  Figure 
c 


7.  Design  Approach  and  Analysis  for  Elliptical 

Airframe 

The  design  technique  and  analysis  for  elliptical 
airframe  is  the  same  as  that,  in  section  II. D. 4.  Again, 
for  purposes  of  analysis,  a  CSMP  program  was  written 
(Appendix  E)  using  equations  of  aerodynamic  model  and 
control  law. 

Figure  2.18  shows  the  yaw  acceleration  response  of 
the  elliptical  airframe.  The  response  has  0.35  second  time 
constant,  overshoot  about  6%  and  steady  state  error  0.01 
which  are  according  to  the  requirerraents ,  (Appendix  C) . 

Figures  2.19  and  2.20  show  the  required  body  yaw 


angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response. 


040  040  0.00  0.70  1.00  140  1.00  1.70  7.90  240  2.00  2.70  3.00  340 

TtSCCI 

Fig.  2.18  Yaw  Acceleration  (ny)  vs  Time  (t) 

Uncoupled  Yaw  channel 

Elliptical  airframe  (lGee  Command.  a„ 


OJO  OJi  OJO  0.70  1.00  lJt  1.00  1.70  2  JO  321  3.K  2.71  2.00  JJI 

T1SCCI 

Fig.  2.19  Yaw  Angular  Rate  (r)  vs  Time  (t) 
Uncoupled  Yaw  Channel 

Elliptical  airframe  (lGee  Command,  a  «0) 
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LINEAR  UNCOUPLED  ROLL  CHANNEL 


Generally,  the  purpose  o£  a  roll  channel  autopilot  is 
to  command  the  ♦  signal,  in  order  to  roll  the  missile 
through  an  angle  ♦  ,  measured  from  the  vertical  axis. 

The  roll  channel  of  a  coordinated  bank-to-turn 
autopilot  is  commanded  to  roll  the  missile  so  as  to  put  the 
preferred  maneuver  direction  of  the  missile  in  the 
direction  of  the  guidance  acceleration  command,  while  the 
pitch  channel  acceleration  is  commanded  to  produce  the 
total  magnitude  of  the  guidance  acceleration  command. 

The  desired  maneuver  plane  acceleration  should  be 
attained  as  rapidly  as  the  achieved  body-fixed  pitch 
acceleration.  To  accomplish  this,  the  uncoupled  roll 
channel  autopilot  was  designed  to  have  the  roll  angle  time 
constant  equal  to  the  time  constant  of  the  normal 
acceleration  achieved  by  the  uncoupled  pitch  channel 
autopilot.  This  is  not  necessarily  the  optimal 
relationship  betweem  the  two  time  constants. 

1.  Roll  Aerodynamic  Model  and  Control  Law 

A  block  diagram  of  the  uncoupled  Roll  channel  is 
shown  in  Figure  2. 21, taken  from  [Ref.  2].  In  this  diagram 
the  aerodynamic  model  is  involved.  Also,  the  roll  control 
law  is  shown  in  Figure  2.22.  These  figures  were  taken  from 
[Ref.  2]. 

The  roll  control  law  is  commanded  by  roll  angle  4>c 

is  governed  by  roll  angular  rate  p  and  roll  angle  $  and 
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determines  the  required  command  5R  to  an  actuator  which  is 

c 

approximated  as  a  first  order  lag  at  30Hz  (188.4  rad/sec) . 
The  aerodynamics,  linearized  about  a  trim  angle-of-attack 
(ag)  equal  to  zero. 

2.  Equations  of  Aerodynamic  Model  for  Circular 
Airframe 

The  aerodynamic  model  is  represented  by  the 
following  equations: 

p-  •Sq  (radlsec1)  di.E.2-1) 

cfc  -  PCradlsec) 


(II. E. 2-2) 


where: 


q  =  1650  (lbs/ft^) 
s  =  3.14159  (ft2) 
d  »  2  (ft) 

Ixx  =  40  (slug-ft2) 


Cl  =  0.031 
6R 


As  referred  to  in  Appendix  A  and  Table  I. 
Substituting  the  above  values  in  equation  II. E. 2-1,  it 
yields: 

P  -  S. 03M 8  •  <5p  (ii. e. 2-3) 

3.  Equations  of  Aerodynamic  Model  for  Elliptical 
Airframe 

The  aerodynamic  model  equations  are  similar  to 
those  for  the  circular  airframe,  with  the  only  differences 
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0.023 


in  the  constant  Ixx=110 (slug-ft  )  and  C( 


Hence  they  are: 

p  -  a. -Sr  CtQLci(se<?')  (ii.E.3-1) 

4>=  P  OradlseO  (ii.E.3-2) 

4.  Equations  of  Control  Law  for  Circular  and 
Elliptical  Airframe 

The  yaw  control  law  for  circular  and  elliptical 
airframe  is  shown  in  Figure  2.22. 

a.  Roll  Angle  Compensator  Equation 

*2.2 


Xr 


9 


+  1 


(II. E. 4-1) 


inverse  Laplace  transforming  the  equation  II. E. 4-1 
it  becomes: 

X  =  (4>c-4>)  (II.  E.  4-2) 

b.  Rate  Compensator  Equation 

0.15065-  (^*+0 

(4+0 

Utilizing  inverse  Laplace  transformation,  and 


v= 


(x-p) 


(II. E. 4-3) 


minor  manipulating  the  equation  II. E. 4-3,  one  obtains: 

V  r -5- Y+ 50. 333XH0*  (X-P)  + 0.1555  (x-p)  (ii.e.4-4) 


c.  Psuedo  Differentiator  Equation 

0.01303 


(II. E. 4-5) 


Using  inverse  Laplace  transformation  and 
rearranging  the  equation  II. E. 4-5,  it  becomes: 

X,  -  -6X4  +  0.0132-  P 


(II. E. 4-6) 


d.  Equation  of  Actuator  Compensator 

This  network  of  roll  control  law  is  different 
for  circular  and  elliptical  airframe. 

(1)  Circular  Airframe 


&♦*(**/ (I—> 

Using,  inverse  Laplace  transformation  and 

rearranging  the  equation  II. E. 4-7,  it  yields: 


4’c  -  -i5jRc  +  0.'l5^3(Y-X<)  +  <5(Y-X>)  (ii.e.4-8) 

(2)  Elliptical  Airframe 

^  (y  x  ■)  UI-E.4-9) 

■nr + * 

Inverse  Laplace  transforming  and 
rearranging  the  above  equation,  it  becomes: 


Sec  r  -isS’Bc  +0.56993  (  Y-3t,  )+6J.5«(y-xO 


(II. E. 4-10) 

e.  Actuator  Equation 

(II. E. 4-11) 

— i-  +  i 

Utilizing  inverse  Laplace  transformation,  the 
equation  II. E. 4-11  becomes: 

-  —  A8S*^I  •  <5* +  1 WF95 -32  ( 1 1 .  E .  4-1 2 ) 

Utilizing  state-space  representation,  all  the 


aforementioned  equation  of  the  aerodynamic  models  and 


control  laws  for  the  circular  and  elliptical  airframes,  can 
be  modelled  in  the  state  variables  systems  given  below: 

(1)  For  the  circular  airframe 

:  o  o  o  o  o  o  g.owl  I  p  |  |  o  I 


o  o 
iu  -8 


Z  -W56  0.?86 

o  o 


o  0  o  o  X  -1U 

-5  0  0.353  -<W0*W  *  Y  +  -0.8W  $c 


p 

• 

X 

• 

Y 

- 

X, 

Sc 

Sr 

0  o  0  0  “6  O  6.2S3  0 

i 

-<M03  0.(11  0.M83  (<1.3(9  -M.«  “5  "0.03?  ^  -0.121 

<D  O  O  O  O  (079133  O 

(II. E. 4-13) 

(2)  For  the  elliptical  airframe 

o  O  o  o  o  o  2.H«  p  O 
iOOOOOO<|>  O 

0  -8  OO  O  O  X  -1H 

-  -0.T56  0.8^6  0  -5  0.353  0  -0.4044  *  Y  +  -0.856  4>c 

O  0  0  0  -60  6.283  X(  0 

-<P4*  -0504  0.702  -3.(41  59.(40  -15  -a 588  -aS04 

0  0  0  O  O  (0T913  <JR  O 


(II. E. 4-14) 


where  the  state  variables  are: 

P  :  roll  angular  rate 
♦  :  roll  angle 

X  :  output  of  roll  angle  compensator 
Y  :  output  of  rate  compensator  network 
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X^  :  output  of  pseudo  differentiator 


6R  :  input  command  in  the  actuator 
c 

sR  :  roll  tail  incidence 

The  state  «,XrPfYfX1#  are  shown  in  Figure 

c 

2.22. 

5.  Design  Approach  and  Analysis  of  Uncoupled  Linear 
Roll  Channel  for  Circular  and  Elliptical  Airframe 
The  design  technique  and  analysis  is  the  same  as 
that  which  was  used  for  the  Pitch  and  Yaw  channel. 

For  purposes  of  analysis,  two  CSMP  program  were 
written  (Appendix  F)  using  the  state  variable  system 
represented  by  the  aerodynamic  model  and  control  law. 

Block  diagram. 

Figures  2.23  shows  the  roll  angle  response  of  both 
elliptical  and  circular  airframe.  The  time  constant  of  the 
roll  angle  response  is  0.55  seconds,  the  overshoot  3%  and 
the  steady  state  roll  angle  error  equal  to  zero. 

Figure  2.24  shows  the  required  body  roll  angular 
rate  for  both  elliptical  and  circular  airframes.  Only  the 
roll  tail  angular  deflection  figures  2.25  and  2.26  is 
different  for  the  airframes,  due  to  the  method  of 
compensating  for  a  reduction  in  the  elliptical  aerodynamic 
roll  gain.  This  is  desirable  for  the  aerodynamic  roll  gain 
(57.3qsdC  /I  )  to  as  ^ar9e  as  possible  to  minimize 


PHI (RAO) 


o.m  o.a  o.ao  o j*  i.or  l.a  t.so  t.a  3.00  ijs  3. so  3.7s  3.00  3 

TiSCC) 

Fig.  2.24  Roll  Angular  Rate  (p)  vs  Time  (t) 
Uncoupled  Roll  Channel 
Circular  or  Elliptical  airframe 
(1  RAD  Command,  a  =0) 


Fig.  2.25  Roll  tail  incidence  (5R)  vs  Time  (t) 


Uncoupled  Roll  Channel 
Circular  airframe  (1  RAD  Command, 
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Fig.  2.26  Roll  tail  incidence  (5„)  vs  Time  (t) 
Uncoupled  Roll  Channel1' 

Elliptical  airframe  (1  RAD  Command,  ag=0) 
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control  surface  motion. 


It  is  also  desirable  to  have  as 


large  as  possible  to  minimize  the  effects  of 

aerodynamic  control  cross  coupling.  The  circular  airframe 
has  a  considerably  larger  aerodynamic  roll  gain  due  to  a 
much  smaller  roll  inertia  and  a  larger  control  derivative 
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III.  LINEAR  DESIGN  AND  ANALYSIS  OF  COUPLED  AUTOPILOTS  FOR 
CIRCULAR  AND  ELLIPTICAL  AIRFRAMES 


A.  GENERAL 


The  linear  design  and  analysis  technique  began  with 
uncoupled  autopilot  channels.  The  uncoupled  autopilot 
design  technique  was  classical,  using  a  combination  of 
frequency  response  and  root  locus  techniques  [Ref.  2],  to 
achieve  practical  bandwidths  (i.e,  sufficient  high 
frequency  attenuation)  and  in  turn  provides  the  range  of 
required  missile  body  angular  rates  and  control  motions,  as 
mentioned  in  section  II. 


In  this  section,  CBTT  control  laws  are  analyzed  to  add 
control  coupling  for  coordinated  missile  motion.  A  measure 


of  sideslip  control  is  obtained  by  applying  a  roll  angle 
command  to  the  linearized  CBTT  autopilot.  The  relative 
stability  of  the  autopilot  branches  and  means  for  improving 
stability  are  discussed.  An  examination  of  the  autopilot 
sensitivity  to  aerodynamic  cross-coupling  is  made. 


B.  AERODYNAMIC  MODELS  FOR  THE  CIRCULAR  AIRFRAME 

The  aerodynamic  models  of  pitch  and  roll-yaw  channels 
for  the  circular  airframe  are  shown  in  Figures  3.1  and  3.2, 
taken  from  [Ref.  2].  The  models  are  linearized  and  were 
developed  for  stability  studies  in  the  frequency  domain. 
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The  assumptions  for  linearization  are  the  same,  as  those  in 
section  II.C.l. 

Pitch  and  Lateral  channels  are  coupled  via  constant 
missile  role  rate  PQ.  Also,  the  linear  model  has  inertial, 
kinematic  and  aerodynamic  cross  couplings,  which  will  be 
discussed  further,  in  the  non-linear  model. 

1.  Pitch  Channel 

The  linear  aerodynamic  model  of  the  pitch 
channel  consisted  of  the  uncoupled  model  mentioned  in 
chapter  II,  which  is  coupled  with  the  Roll-Yaw  channels  via 
constant  roll  rate  { Pe ) ,  as  shown  in  Figure  3.1. 

In  this  work,  the  Pe  was  set  to  zero,  because 
when  Pe=0,  the  Qe  (i.e  equilibrium  pitch  rate)  has  been 
found  to  have  negligible  influence  in  the  lateral  model 
(i.e  roll-yaw) . 

The  equations,  which  represent  the  pitch 
aerodynamic  model  given  in  [Ref.  2]  are  the  same  as  the 
equations  of  the  uncoupled  linear  model,  mentioned  in 
chapter  II. 

2.  Lateral  Channels 

The  linear  lateral  aerodynamic  model  consisted 


of  the  uncoupled  roll  and  yaw  aerodynamic  channel,  coupled 


via  the  linearized  aerodynamic  derivatives  C. 


Cn  .  Also,  it  is  coupled  with  the  pitch  dynamic  model  via 
6r 
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Linear  Pitch  channel  dynamic  model 


the  Pe,  as  previously  mentioned.  The  lateral  dynamic  model 
is  shown  in  Figure  3.2,  as  given  in  [Ref.  2], 

To  find  the  equations  which  represent  the 
lateral  dynamic  model,  Pe  and  Qe  are  set  to  zero  because  of 
their  negligible  influence  in  the  lateral  channel. 

The  aforementioned  equations  are  as  follows: 


Ayr  q-£  £(  /  W 

(II. B. 2-1) 

l\  Z.  6*  Cy$  4-  Cyjy.  £y  (<Jeg) 

(III. B. 2-2) 

where : 

t  r  57.3 q -Ej  /  Iii  (deglset2) 

(III. B. 2-3) 

E}  r  &  C i\j?  +  .  £r  4  CriSy  •  Sy  (dey) 

(III. B. 2-4) 

P  -  513  3  ^  in  /Ixx  (deg (sec2) 

(III. B. 2-5) 

4>r  P/51.3  (fctdisec) 

(III. B. 2-6) 

h-  (cjy  •  3V  +  Cde§) 

(III. B. 2-7) 

Eqc  Ejf  Ce$.£  (deg) 

(III. B. 2-8) 

#3  f>A/5f.3-r+V<-(qiEilw)  (dtglsct) 

(III. B. 2-9) 

where : 

k 

=  0.48 

S 

«  "(ft)2 

d 

=  2. 0  ( f t) 

q 

*  1650  ( lbf/ft2) 

w 

=  2525  (lbf ) 

Ixx 

=  40.0  (slug-ft2) 

Izz 

=  810.0  (slug-ft2) 

V 

=  0.22 
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=  constant  pitch  rate  (deg/seci 

,2  Linear  lateral  (Roll/Yaw)  dynamic  model 


"3 


-0.019 


n«y  «  -0.053 
Ci 

6y  =  -0.016 
CA  . 

3  =  -0.009 

«R  =  0.035 
A  =  10° 

As  referred  to  in  Appendix  A  and  Table  I. 

Substituting  the  values  into  the  above  equations, 
they  yield: 

T  :  -13.93M-6-M3.W-Jp-  3S.Sf-Sy (d€5ljetzJ(iii.B.2-i0) 

p  3  -133.66 -6+519. IM-St  -23T618  iy(tyUt)  (HI  2-ID 

8  ZO.OU M5  P-T-0.080S  G+O.lil  SyCdez/sec)  (Hi-B*2-12) 

<*>-  P/5T-3  (VCCdlseO  (III. B. 2-13) 

C.  CBTT  AUTOPILOT  CONTROL  LAWS  FOR  CIRCULAR  AIRFRAME 
The  control  laws,  which  were  used  by  the  CBTT 
autopilots  of  the  circular  and  elliptical  airframe,  are 
shown  in  Figures  2.6  and  3.3,  as  given  in  [Ref.  2], 

The  pitch  control  law  for  the  circular  airframe  is  the 
same  as  determined  in  the  uncoupled  pitch  channel  study 


CBTT  lateral  control  law 


section  II. B. 2,  and  therefore,  the  equations  which 
represent  it,  are  the  same  as  in  section  II. B. 3.,  except 
the  equation  for  pitch  tail  incidence  angle  which  is: 

cTp=  -18S.4 -5p 440495.3 -<5>c  (demise c)  (iu.c.i) 

The  uncoupled  roll  and  yaw  laws  are  coupled  via  a 
cross-coupling  branch,  shown  in  solid  line  in  Figure  3.3, 
in  order  to  consist  the  linear  coupled  lateral  control  law. 

In  the  roll  actuator  of  the  lateral  control  law,  the 
frequency  has  been  changed  from  110  rad/sec  to  60  rad/sec, 
for  improving  of  autopilot  stability. 

The  autopilot  cross-coupling  branch  between  Roll  and 
Yaw  channels  has  been  added  to  provide  coordinated  motion. 
Coordinated  motion  or  zero  sideslip  angle  is  achieved  by 
directing  the  body  fixed  pitch  axis  of  the  missile  of  the 
missile  velocity  vector,  so  that  there  is  no  component  of 
missile  velocity  along  the  body  fixed  yaw  axis  of  the 
missile.  Figure  3.4,  as  given  in  [Ref.  2]  shows  the 
attitude  of  the  missile  body  with  respect  to  its  velocity 
vector  V. 

When  commanding  an  upward  maneuver  (i.e,  *  =0 ) ,  the 
missile  body  moves  upward  with  its  pitch  axis  directed  at 
the  velocity  vector  until  it  reaches  the  desired  maneuver 
level  or  angle-of-attack.  No  roll  motion  is  required  to 


maintain  coordination  for  this  maneuver. 
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For  maneuvers  in  the  *=45  or  90  degree  directions,  the 
missile  body  rolls  the  desired  ♦  and  moves  its  pitch  axis 
directed  at  the  velocity  vector,  until  it  achieves  the 
desired  maneuver  level.  For  maneuvers  with  *  higher  than 
90°,  the  missile  must  roll  about  the  velocity  vector,  while 
the  Yaw  channel  directs  the  pitch  axis  towards  the  velocity 
vector  for  minimum  sideslip  angle. 

The  missile  is  forced  to  roll  about  its  velocity 
vector,  when  the  desired  maneuver  direction  is  in  the 
negative  angle-of-attack  direction,  because  it  is  desired 
to  avoid  negative  angles-of-attack.  These  maneuvers  are 
illustrated  in  Figure  3.4. 

The  coordinating  command  rc,  is  a  yaw  angular  rate  from 
a  rate  gyro.  The  command  is  equal  to  -  Pae,  where  ae=a 
(Figure  3.3)  in  the  coordinating  branch  is  the  equivalent 
angle-of-attack,  which  is  exactly  equal  to  angle  of  attack 
(a) ,  in  the  linear  studies.  The  equations  which  represent 
the  lateral  channel  are  as  follows: 


Inverse  Laplace  transforming  and  rearranging 


the  equation  (III. C. 1-2),  it  becomes: 


Y,  -  -5  V  0.0  503  (X-P/5T-3)  +  0.TS5  (X-  P/5T-3)  (III. C. 1-3) 


c.  Psuedo  Differentiator  Equation 


(  Pi 51.5) 


(III. C. 1-4) 


Xt- - Q^iLiS - (  p,„3)  (III.C.1- 

-lr  +  1 

Utilizing  inverse  Laplace  transformation  and 


minor  manipulating  the  above  equation,  one  has: 

X<  r  ~<o\  +  0.073  (  Pl  515) 

d.  Actuator's  Compensator  Equation 


(III. C. 1-5) 


The  frequency  of  the  compensator  was  changed 
from  110  rad/sec  to  60  rad/sec  for  purposes  of  stability, 
as  referred  in  section  III.C.  The  equation  represents  the 
network  is: 


+  i 


(III.C. 1-6) 


Using  inverse  Laplace  transformation,  the 


equation  III.C. 1-6,  one  obtains: 


f.  =  -i5  Src  +  0.25  (Y-X,)  +15  (y-Xi)  (ui.c.i-7) 

e.  Actuator  Equation 


—  c5rc 

mA  *  * 


(III.C. 1-8) 


Utilizing  inverse  Laplace  transformation  and 
rearranging  the  equation  III.C. 1-8,  it  yields. 


S.  z-  itt.  4  S*  +  WS5.  S2  Arc 


p.  --1»#-TUKT  ,w  I  (III. c.  1-9) 

All  the  referred  to  above  state  variables,  are 
shown  in  figure  3.3. 

2.  Yaw  Channel 


a.  Acceleration  Compensator  Equation 


y  _  Q.  32 

<m+  i 


tly 


(III. C. 2-1) 


Inverse  Laplace  transforming,  equation 
I I I. C. 2-1  becomes: 

V  -  -5Y41.6  t|y  (hi. c. 2-2) 

b.  Actuator  Compensator  Equation 

5ycr  (y-0.M58a.PfeT3+t)  (III. C. 2-3) 

Utilizing  inverse  Laplace  transformation  and 

rearranging  the  equation  III. C. 2-3,  it  yields: 

•  # 

QYC  =O.lW5(Y+Yk?.3-0.«S  a-P/5T.3)4  M-85(Y+r|5T.3-a«S-a*P/5T.3) 

(III. C. 2-4) 

c.  Actuator  Equation 

<5y=-l$*.H-$Y  +  40?95.3  •  5yc  (III. C. 2-5) 

Making  use  the  state-variable  representation. 


the  above  equations,  of  the  aerodynamic  models  and  control 
laws,  can  be  modelled  in  the  following,  state  variable 
equation  system: 


O  O  ©  -13.W  o  o 

0  0  0  -15M6  0  0 

aoiw  o  o  o  o  o 

-1  0.W5  o  -aow  o  o 

o  0  “JW  0-8  0 

O  OOfJa  *0.116  ailTV  03535 
o  o  o  -o.m  o  o 

o  00053  -0225  0.WW  0.0«J»  J5.T5 
ooo  o  0  O 

o  o  o  -azo  o  o 

o.o<w*ooo«  O  -0.1W  0  0 

o  o  o  o  o  o 


O  HMH  O  O  -36-*?  V 

0  529-19  o  0  p 

o  o  o  o  o  4> 

o  0  0  o  0JK8  € 

0  0  0  O  0  * 


o  o  -aw  o  o  o.iott  ,  ^ 

-6  0  O.?0M  0  O  -0.315  *4 

-13.6  -15  0.131  0  °  Ak 

0  4H913  HIM  ooo  fa 

0  0  0  *5  0  0  V 

o  O  o  -a  482  3.31S  MW  ** 

oooo  $y 


0 

iu 

0.8JC  (|>c 
0 

0.2115 


(III. C. 2-6) 


where  the  state  variables  are: 


& 


yaw  angular  rate 
roll  angular  rate 
roll  angle 
sideslip  angle 

output  of  roll  command  compensator 
output  of  roll  rate  compensator 
output  of  pseudo  differentiator 
input  command  in  roll  actuator 


roll  tail  incidence 

output  of  yaw  acceleration  compensator 
input  command  in  yaw  actuator 
yaw  tail  incidence  angle 


All  the  referred  before  state  variables  are 


shown  in  Figures  3.2  and  3.3. 


OW5*3 
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D.  ANALYSIS  OF  LINEAR  CBTT  AUTOPILOTS  FOR  CIRCULAR 

AIRFRAME 

In  the  analysis  of  the  linear  CBTT  autopilots  it  is 
assumed  that  the  missile  is  initially  in  the  desired 
maneuver  plane,  is  trimmed  at  ten  degrees  angles-of-attack, 
the  equilibrium  roll  rate  Pe  and  equilibrium  pitch  rate  Qe 
were  set  equal  to  zero. 

For  purposes  of  analysis,  a  CSMP  program  was  written 
(Appendix  G)  using  the  equations  of  the  aerodynamic  models 
and  control  laws  of  the  linear  CBTT  autopilots  for  the 
circular  airframe. 

Figure  3.5  shows  the  sideslip  angle  for  the  linearized 
CBTT  autopilot,  having  a  maximum  value  of  2.32°.  Since  the 
model  is  linear,  and  therefore  the  magnitude  of  the 
sideslip  is  directly  proportional  to  the  magnitude  of  the 
input  command,  a  one  radian  roll  angle  command  was  used  for 
convenience. 

The  coupled  autopilots  have  a  roll  angle  response  which 
is  the  same,  as  the  uncoupled  roll  channel  response, 
differing  slightly  in  the  overshoot,  as  shown  in  Figure 
3.6. 

The  aforementioned  values  and  plots  are  exactly  the 
same  as  the  plots  given  in  the  [Ref.  2]  therefore  its 
results  are  verified. 
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Pig.  3.5  Sideslip  angle  (g)  vs  Time  (t) 

Linear  CBTT  Autopilot  (1  RAD  Roll 
Command);  Circular  airframe  (a  »10  deg, 
Pe»Qe»0  deg)  e 
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Comparing  the  state  variable  form  system  with  these 
previous  discussed  in  the  section  II. D  and  II. E,  this  is 
the  twelfth  order  system  because  it  includes  roll  and  yaw 
autopilots  and  the  matrices  A,B  are  different  than  the  A,B 
matrices  at  the  uncoupled  roll  and  yaw  autopilot,  because 
of  the  existance  of  the  coordination  branch. 

E.  AERODYNAMIC  MODELS  FOR  ELLIPTICAL  AIRFRAME 

The  dynamic  models  for  elliptical  airframe  are  the  same 
as  those  of  circular  airframe  and  are  shown  in  Figures  3.1 
and  3.2. 

The  differences  of  the  circular  and  elliptical  models 
are  in  the  following  constants  and  linearized  aerodynamics 
derivatives  s 

«  2475  (lbf) 

*  110  (slug-ft2) 

*  790  (slug-ft2) 

-  0.015 

i 

-  -0.054 

=  0.014 

-  0.024 

*  -0.039 


w 

Ixx 

Izz 
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*  -0.010 


Ct 

6  a  -0.027 

C % 

«R  =  0.023 


as  referred  to  in  Appendix  A  and  Table  I.  The  equation 
which  represent  the  aerodynamic  model ,  are  as  follows: 


Y  Z  +  -fo- 39.3242. Jv 

p  r  -5M.04-5v 

$  -  P/S7-3 

^  z  0.1*452  p- T  +  0.0643W  4-0.0150f'Sy 


(III. G. 1) 
(III.E.2) 
(III. E. 3) 
(III.E.4) 


F.  CBTT  AUTOPILOTS  CONTROL  LAWS  FOR  ELLIPTICAL  AIRFRAME 
The  control  laws,  which  were  used  by  the  CBTT 
autopilots  of  the  elliptical  airframe  are  the  same  as  those 
of  the  circular  airframe  and  are  shown  in  Figure  2.6  and 
3.3. 

The  differences  which  exist  in  the  control  laws  of  the 
two  airframes  are  the  gain  kYP  of  the  coordinated  branch, 
and  the  gains  of  the  actuators. 

Hence,  the  equations  that  represent  the  control  laws 

are: 

1.  Pitch  Channel 

The  linear  pitch  control  law  equations  for 
the  elliptical  airframe,  are  the  equations  of  the  Pitch 

uncoupled  channel,  as  given  in  section  II. C. 5. 
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2.  Roil  Channel 


Same  equations  as  in  par.  III.C. 


a. -x  =-SX+i7.6  (II. C. 1-1) 

b.  Y<r -5 Y,+ 0.0505  (X-P|ST.5)^0.?55 (X-P/5T.5)  (II. C.  1-3) 

c.  -OX.,  +  0.075  (  P/57-3)  (II. C. 3-5) 


d.  Actuator  Compensator 

( V: 


(II. C. 3-5) 


(III. P.2-1) 


Utilizing  inverse  Laplace  transformation  and 
rearranging  the  equation  III. F. 2-1,  it  becomes: 

=  -15^  + 10^5 (VX|) 462.53  (Y,-Y,)  (m.p. 2-2) 

e.  £*=  -138.4  £*  +  J.0.?055XIo3Src  (III.C. 3-8) 


3.  Yaw  Channel 

a.  Acceleration  Compensator  Equation 


yj  0.839  (III. F. 3-1) 

V._„ 

Inverse  Laplace  transforming  and  rearranging 


the  equation  III. F. 3-1,  it  yields: 

Y  =  -4Y  +  3.350  Hy 


(III. F. 3-2) 


b.  Actuator  Compensator  Equation 

Jyc  -  ^0JSO9±±±O9  (  y  +  y/5T. 3 . 5.  P/5f  3) 

5 

(III .F. 3-3) 

Utiliing  inverse  Laplace  transformation  and 
rearranging  the  above  equation,  one  obtains: 

$Yt  =  «•«*(>+ £5*  f&)  +4-OI  (Y*S}-55  )  (XII.F.3-4) 


10' 


c.  Actuator  Equation 

Same  as  in  par.  III.F.2.e 

Utilizing  the  state-variable  representation, 
the  equation  which  represent  the  aerodynamic  model  and  the 
control  laws,  can  be  modelled  in  the  following 
twelfth-order  system: 
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(III. F. 3-5) 


?re 


state  variables  are: 
yaw  angular  rate 
roll  angular  rate 
roll  angle 
side  slip  angle 

output  of  roll  command  compensator 
output  of  roll  rate  compensator 
output  or  psuedo-dif f erentiator  network 
input  command  in  roll  actuator 


3 


R 


roll  tail  incidence  angle 


¥ 


output  of  acceleration  compensator  in  yaw 
control  law 


G.  ANALYSIS  OF  LINEAR  CBTT  AUTOPILOTS  FOR  ELLIPTICAL 

AIRFRAME 

In  the  analysis  of  the  linear  CBTT  for  elliptical 
airframe,  the  same  assumptions,  which  were  made  for  the 
circular  airframe,  are  valid. 

A  CSMP  program  was  written  (Appendix  G)  with  the  above 
equations,  which  represent  the  aerodynamic  models  and  the 
control  laws  of  the  linear  CBTT  autopilots. 

Figure  3.7  shows  the  sideslip  angle  S  for  the 
linearized  CBTT  autopilot,  having  maximum  values  of  -0.46° 
and  0.26°. 

Figure  3.8  shows  the  roll  angle  response  which  has  an 
overshoot  of  4  percent  and  is  about  the  same  with  the  roll 
angle  of  uncoupled  roll  autopilot.  The  above  plots  verify 

the  results,  given  in  the  [Ref.  2]. 
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Fig.  3.7  Sideslip  angle  (8)  Ys  Time  (t) 

Linear  CBTT  Autopilot  (1  RAD  Roll 
command)  Elliptical  airframe  (a  *10°, 
Pe=Q  =0°) 


H.  RELATIVE  STABILITY  OF  CBTT  AUTOPILOTS  OF  CIRCULAR  AND 


ELLIPTICAL  AIRFRAMES 


When  the  coordinating  branch  gain  Kyp  was  set  to  zero 
in  CBTT  autopilot  of  the  circular  airframe,  the  maximum 
sideslip  angle  increased  to  2.8  degrees.  Therefore,  for 
the  circular  airframe  the  coordinating  branch  is  not  very 
effective  in  helping  the  yaw  autopilot  to  reduce  sideslip 
angle . 


If  k^p  is  increased  to  unity,  as  it  is  for  the  CBTT 
autopilot  of  the  elliptical  airframe,  the  maximum  sideslip 
remains  below  one  degree  for  3  seconds  but  the  autopilot  is 
unstable.  Modifying  of  the  circular  airframe  will  be  done 
in  the  non-linear  CBTT,  for  minimization  of  the  sideslip 
ang le . 


Removal  of  the  aerodynamic  cross-coupling  (i.e  C 


'6 


C^  ,  Cn  )  had  no  effect  on  the  sideslip  angle.  This 
V  *R 


showed  that  the  aerodynamic  cross-coupling  plays  an 
indirect  role  in  determining  the  quality  of  sideslip 
control  by  determining  the  relative  stability  of  the 
coordination  branch  or  the  magnitude  of  the  coordination 
gain  K^p.  The  magnitude  of  sideslip  angle  is  dependent  on 
the  nulling  effects  of  two  parallel  paths  shown  in  Figure 
3.9,  as  given  in  [Ref.  2]. 


The  contribution  of  yaw  acceleration  n^  to  the  maximum 
0  is  negligible  and  therefore  neglected.  6  is  formed 
mainly  by  the  subtraction  of  the  kinematic  paths  of  a@p  and 
the  yaw  angular  rate  r. 

The  coordination  is  obtained  in  the  CBTT  control  law  by 
commanding  the  yaw  autopilot  with  a  yaw  angular  rate 
command  rc  (Figure  3.2)  of  Kyp  a0p ,  which  forces  r  to  be 
equal  to  ae.p  and  therefore  nulling  as  shown  in  Figure 
3.9.  The  nulling  process  will  be  accomplished  more 
efficiently  if  Kyp  =1.0.  The  reason  the  sideslip  is  not 
nulled  completely  is  that  ny  is  not  zero  and  r  can  not 
equal  r  over  all  frequencies. 

The  linearized  CBTT  autopilot  of  the  circular  airframe 
was  unstable,  when  the  uncoupled  channel  control  laws  used 
with  a  coordinating  branch  gain  Kyp  of  unity.  The 
autopilot  was  stabilized  by  decreasing  to  0.458.  The 
actuator  of  roll  control  channel  has  a  frequency  of  60 
rad/ sec. 

The  stability  of  the  CBTT  autopilot  for  the  elliptical 
airframe  is  considerably  better,  due  to  the  setting  of  Kyp 
at  unity,  which  results  in  a  substantial  decrease  in  the 
magnitude  of  sideslip  angle  as  shown  in  Figure  3.9. 
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I.  CONCLUSIONS 


1.  The  linear  coupled  autopilots  verify  the  results  of 
the  uncoupled  autopilots  and  keep  the  sideslip  angle  small 
for  coordination  motion,  (section  II. D.  and  II. E) 

2.  Further  Conclusions  about  the  relative  stability  of 
the  autopilots,  about  the  gain  Kyp  of  the  coordination 
branch  and  about  the  sensitivity  of  the  autopilots  to  the 
aerodynamic  cross-coupling  are  given  in  chapter  7  of  [Ref. 


IV.  NONLINEAR  ANALYSIS  OF  CBTT  AUTOPILOT  FOR  CIRCULAR  AND 
ELLIPTICAL  AIRFRAME 


A.  GENERAL 

The  three  dimensional  nonlinear  aerodynamic  models 
used  for  the  analysis  are  shown  in  the  Figures  4.1  and  4.2. 
For  the  configuration  under  consideration,  the  assumptions, 
which  have  been  found  to  be  consistent  are  referred  to  in 
chapter  5  of  [Ref.  2]. 

Nonlinear  aerodynamics  and  mass  parameters  values  are 
presented  in  Appendixes  A  and  B,  as  taken  from  the  [Ref. 

2].  The  same  flight  conditions  (i.e  60kft  altitude  and 
M=»3. 95)  as  in  linear  CBTT. 

In  this  section  the  control  laws  for  the  non-linear 
CBTT  autopilots  are  discussed.  The  results  for  commanding 
the  CBTT  autopilots  for  both  airframes  are  shown  and  the 
desired  aerodynamic  model  to  enhance  CBTT  performance  is 
determined.  The  performance  of  the  CBTT  autopilot  for  the 
elliptical  airframe  is  determined,  when  no  lateral 
aerodynamic  cross-coupling  exists  and  with  ideal  airframe 
dynamics.  Also  the  performance  of  the  ?qCBTT  autopilot  for 
the  circular  airframe  is  determined,  when  no  inertial  and 
kinematic  cross-coupling  exist  into  the  pitch  autopilot. 

The  conclusion  of  the  nonlinear  analysis  are  stated. 
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B.  CONTROL  LAWS  FOR  ELLIPTICAL  AIRFRAME 

The  control  laws  used  for  the  following  nonlinear  3-D 
studies  were  the  same  as  those  used  for  the  linear  studies 
in  section  III,  which  are  shown  in  Figures  4.3  and  4.4 
taken  from  [Ref.  2] ,  except  for  the  gain  as  shown  in  the 
bold  line  of  the  coordination  branch  in  Figure  4.4.  The 
new  gain  a  is  held  constant  at  one  degree  magnitude  for 
angles  of  attack  less  than  one  degree  positive  and  greater 
than  negative  5  degrees.  For  angles-of-attack  greater  than 
one  degree  positive,  the  gain  a  is  equal  to  the  angle  of 
attack.  This  maintains  coordination  for  very  small 
angles-of-attack. 

Gravity  effects  were  not  included  in  the  linear 
studies,  because  it  was  assumed  to  have  a  negligible 
influence  on  autopilot  stability  and  response  perturbations 
about  a  missile  trim  conditon. 

Gravity  effects  were  included  in  the  nonlinear 
studies,  where  the  missile  body-fixed  yaw  axis  will  be 
subjected  to  the  full  force  of  gravity  and  may  therefore 
have  a  significant  influence  on  sideslip. 

1.  Pitch  Control  Law 

The  pitch  control  law  for  the  elliptical  airframe 
is  the  same  that  is  used  in  the  linear  studies  in  section 
III.F.,  except  in  the  command  applied  to  the  channel. 


CBTT  Nonlinear  lateral  control 


In  inertial  rectangular  coordinates,  the 


I 

* 

I 


I 


acceleration  commands,  are: 

a.  nc  (acceleration  command  in  inertial  Zy 

direction  as  in  figure  2.4)  =  nz  -  cos*  IV. B. 1-1 

where:  nz  :  guidance  command  (gees) 

c 

-cos*  :  anti-gravity  bias  command  (gees) 

0  :  Elevation  Euler  angle  =  (qcos*-rsin*) dt 


b.  ny  acceleration  command  in  inertial  Yv 
c 

direction,  which  is  equal  to  zero,  because  it  is  desired 
coordinated  motion,  in  the  polar  coordinates,  the 
acceleration  commands  are: 


a.  'Vic  r  -  [  ('’Uc-  c«d)*+‘nYc2  ] 

z  -  [  (nu-cos©)2]^ 

-  -  (ti2c-  cose) 


(IV. B. 1-2) 


(IV. B. 1-3 ) 


hence  »c  ■  ni  and  n  *0,1,2... 

Since  the  pitch  control  law  of  the  elliptical 
airframe  does  not  have  an  integration  in  the  acceleration 
error  path  (i.e  in  the  acceleration  compensator  as  in  the 
case  of  circular  airframe)  and  requires  a  gain  in  series 
with  the  acceleration  command  shown  in  Figures  4.3  it  was 
necessary  to  modify  the  antigravity  command  as  follows  to 
assure  an  anti-gravity  bias  of  just  one  gee. 
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nc  (elliptical  airframe  inertial  acceleration  command)  = 

■  nz  -  0.913  cos©  (III.B.1.4) 

c 

The  autopilot  was  tested,  using  a  commanded  2  gee 
(0°,  180°),  that  denotes  an  inertial  command  of  2  gees, 
which  is  first  applied  in  the  0°  or  upward  direction  at  t  = 
2  seconds  for  3  seconds  and  at  the  fifth  second  is  applied 
in  the  180°  or  downward  direction.  Since  at  first  both  the 
missile  roll  angle  and  roll  angle  command  are  at  zero 
degrees,  there  is  no  roll  motion  and  the  missile  turns 
upward  as  a  skid-to-turn  controlled  missile.  At  fifth 
second  the  missile  is  commanded  to  roll  through  180  degrees 
while  moving  in  a  coordinated  manner  in  yaw  and  roll  to 
minimize  sideslip  angle  and  prevent  or  minimize  negative 
angle-of-attack. 

The  other  equations,  which  represent  the  nonlinear 
pitch  control  law  are  the  same  as  for  linear  studies 
section  III.C.  and  the  state  variables  are  shown  in  Figure 

The  equations  are: 

X  r  -150X  +  150  T)*b  (iv.b.1-1) 

*  r  -6Y- 0.5305  tlc  +  0.^sx  (iv.B.l-2) 

Jpc  -  -0.383?  (Y+q|5T.3)-3.0?fY-ql5?.3)  (iv.b.i-3) 

Sf  =  -  i88*M -5p  +  *iOW.3*£f>c  (iv. b.  1-4 ) 


4.3 


where : 


n„  :  achieved  body-fixed  acceleration  in  Z  axis 


nc  :  inertial  acceleration  command. 


2.  Lateral  Control  Laws 

The  lateral  control  laws  are  the  same  as  those 
used  for  the  linear  studies,  except  the  modification  in  the 
coordination  branch. 

The  command  a  for  the  roll  control  law  is  zero 
c 

for  the  first  5  seconds  and  180°  affterwards,  as  derived  in 


IV. B. 1. 


The  equations  which  represent  the  lateral  control 


laws  are  the  same  those  of  section  III.F. 


Ya  r  -q  y2  +-3.356  *WS  (iv.b.2-1) 

=  »“*  "  A-  Sfe) 


{IV. B. 2-2) 


where: 


9*=®  \ 

a,=  a  J 


a>r 


a<l° 

at  =  o  i 

1,  =  -  8Jti  -  «•€  <t>+  iT.t  4>c  (IV.B.2-3) 

Si  z -18 8.4  Si  +  10395.3  Si c  (iv. B. 2-4) 

%  =-5Kf+0.05033Ci,-P|s?.3)  +  MS5(X,-p/5i.j)  (IV.B.2.5) 

X2— ~6Xi  + 0*03?  ( P /5T*3^  (iv. b. 2-6) 

5|>cz  -15 +  1.0425  (Y(-X2l+62-55  Ofl-Xi)  (lv.B.2-7) 


£*  =  4  m<35.3  -Sz C  (IV.B.2-8) 

The  state  variables  are  shown  in  Figure  4.4.  For 

purposes  of  analysis,  a  CSMP  program  was  written  (Appendix 
H)  . 

C.  AERODYNAMIC  MODELS  FOR  ELLIPTICAL  AIRFRAME 

The  aerodynamic  models  of  Pitch  and  Roll-Yaw  channels 
for  the  elliptical  airframe  are  shown  in  Figures  4. 1,4. 2. 
The  airframe  dynamics  under  consideration,  have  a  level  of 
complexity  sufficient  to  determine  the  critical  areas  of 
concern  regarding  the  stability  and  response  of  CBTT 
control . 

1.  Pitch  Channel 

The  nonlinear  aerodynamic  model  of  the  Pitch 
channel  consisted  of  the  uncoupled  model,  coupled  with  the 
Roll-Yaw  channels  via  the  inertial  cross-coupling  of  rp 
into  q  and  the  kinematic  cross-coupling  of  -s.P  into  a. 

Also  the  antigravity  command  cos*  cose  is  inserted  into  nz. 
This  model  is  shown  in  Figure  4.1. 

The  above  terms  consist  of  the  nonlinear  terms  of 
the  pitch  channel. 

The  equations  which  represent  the  pitch 
aerodynamic  model  are: 

z  - q  ^  CN  /  W  (gees) 
q  =  P’f/5?.3  +  51.3  C4  llyy 
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(IV. C. 1-1) 
(IV. C. 1-2) 


(IV. C. 1-3) 
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equations  that  represent  the  aforementioned  model  are  as 


follows : 


“V  - 
6  = 
r  - 

p  = 

<t>  r 


9*3  +  Cvfr  •Si')  1  w 

(IV. C. 2-1) 

Kt»YB  +  a.p|5?.3  -r  (devisee) 

(IV. C. 2-2) 

-  9  *  P  /  5^-3  +  51.3  9  C 

•5*)  (l?^ 

(IV. C. 2-3) 

51.3  9  ^  C  ^^4  *  ^  4  /  Ixx 

(IV. C. 2-4) 

P/5T-3  (-faiUec) 

(IV.C.2-4a) 

P  (deglsec) 

(IV. C. 2-5) 

where: 

K  »  0.48 

q  -  1650  (lbf/ft2) 

S  ■  it  (ft2) 

d  »  2  (ft) 

w  -  2475  ( lbf ) 

Ixx  »  110  (slug-ft2 

I  ,  ■  853  (slug-ft2) 
zz 


as  referred  to  in  Appendix  A. 


Functions  Cy  (a) ,  Cy  (a) ,  Cn  (a) ,  Cn  (a) , 

8y  6  8  4y  ®y 


C.  (a) ,  C.  (a)  and  C_  (a)  are  nonlinear  functions  which 
»  SR  !R 

vary  with  a  angle-of-attack  and  referred  in  Appendix  B. 

Substituting  the  values  into  the  referred  before 


equations  one  obtains: 


1.3333  (CV^  +  Cyjy . 5y)  (iv.c.2-5) 

6  =  0.63998  (crf8  +  Cy^-Jy) +0.0 W52 a- P-r  (iv.c.2-6) 
f  r-0.0W59*P+696.M2fcig  (lv.c.2-7) 

p  -5M00.391?  (C^.g+C^.Jy  +  C^.&O  (iv.c.2-8) 

<£j  =  P  (IV. C.  2-9 ) 

D.  ANALYSIS  OF  NONLINEAR  CBTT. AUTOPILOTS  FOR  ELLIPTICAL 
AIRFRAME 

For  purposes  of  analysis,  a  CSMP  program  was  written 
(Appendix  H)  with  the  aforementioned  equations  of  the 
aerodynamics  models  and  control  laws  of  the  nonlinear  CBTT 
autopilots. 

To  decrease  computer  time,  appropriate  initial 
conditions  were  inserted.  The  initial  conditions  are: 
a  :  angles  of  attack  3  2.41  degrees 
6p  :  pitch  tail  incidence  *  0.658  degrees 
e  :  pitch  Euler  angle  3  0.636  rad 

X  :  output  of  pitch  acceleration  feedback  lag  3  lgee 

Y  :  output  of  pitch  acceleration  compensator  3  -0.0105 

6  :  pitch  actuator  command  3  0.0148  rad 

*c 

The  achieved  maneuver  plane  acceleration  n  is 

z 

calculated  from  the  body-fixed  accelerations  n,  (Figure 

Z0 

4.6)  and  (Figure  4.7)  and  the  roll  angle  4  (Figure  4.13) 
as  follows: 


(IV.D-1) 


Yl?  =  co&4>  +  ityg*  s'na<t» 

Figure  4.5  through  4.16  show  the  responses  of  n_, 

z 

nv  /  a'  Q/  P/  t-/  •/  $  ,  iv  and  fiR.  Also  these  figures 
8  p  y 

and  table  II  show  the  time  constants  and  the  overshoots  of 
the  referred  before  variables. 

Figure  4.5  shows  that  the  achieved  inertial 

acceleration  n„  doesn't  satisfy  the  requirement  of  the 

z 

overshoot  (i.e  below  10  percent)  and  it  is  not  acceptable. 
The  results  of  this  analysis  verify  the  result  of  the  [Ref. 
2). 

E.  CBTT  PERFORMANCE  WITH  IDEAL  AIRFRAME  DYNAMICS  FOR 
ELLIPTICAL  AIRFRAME 

The  purpose  of  the  following  simplifications  to  the 
airframe  dynamics  model  is  to  isolate  the  critical 
cross-coupling  paths,  which  have  caused,  the  transients  in 
the  maneuver  plane  acceleration  responses. 

The  dynamic  model  without  the  coupling  paths  will  be 
referred  to  as  ideal  dynamics.  Although  ideal  dynamics  are 
not  physically  attainable,  it  is  a  useful  goal  for  both 
autopilot  and  airframe  designers. 

The  same  guidance  commands  are  applied  to  the  CBTT 
autopilot  of  the  elliptical  airframe  as  in  par.  IV.B.l 
(i.e,  2  gees  (0°,  180°)  but  with  lateral  aerodynamic 
cross-coupling  removed.  That  means,  the  aerodynamic 


Fig.  4.16  Roll  tail  incidence  (6R)  vs  Time  (t) 

CBTT  of  Elliptical  airframe;  2  Gees  (0 


cross-coupling  was  removed.  The  major  influence  of  the 
aerodynamic  cross-coupling  in  the  elliptical  airframe 
response  has  been  to  decrease  the  overshoot  in  the  achieved 
maneuver  plane  acceleration  resulting  from  the  second 
guidance  command. 

In  addition  to  aerodynamic  cross-coupling  the 
kinematic  cross-coupling  of  -6p  into  a  and  inertial 
cross— coupl ing  of  pr  into  q  were  removed.  Therefore,  the 
only  cross-coupling  which  exist  in  the  airframe  dynamic 
model  are  the  kinematic  coupling  of  ap  into  8  and  the 
inertial  coupling  of  -qp  into  r.  There  is  also  the 
autopilot  cross-coupling  of  the  coordinating  command  from 
the  Roll  to  Yaw  channel. 

The  equations,  which  represent  the  modified 
aerodynamic  model  by  removing  the  referred  cross-coupling 
are  same  with  those  as  in  paragraphs  IV.B.l,  IV. B. 2,  IV.C.l 
and  IV. C. 2,  with  the  only  difference  in  the  following 


equations: 


1. 


and 


Pitch  Aerodynamic  Model 

^  Z  5T.3  9  Q*  ttyv  C  <U.g  /sec  ) 
Z  *51. *541  CM 


&  -  s+ 

i  c)  -  (  k  2j  £  Cm  Iw) 

Z  q  -  0.639*8  Qv  (demise c) 


(IV.  E.  1-1 ) 


(IV. E. 1-2) 


2.  Lateral  Aerodynamic  Model 


V  =  -~~-4  (Cyi* .6+0nfy.5y) 

“♦•3  422 

-  -a0W52  «|«p  +696.m?  (<*r8+C*Sf-S'i)  (^Jsec2) 

(IV. E. 1-3) 


•  5i.3  q  & 

V  "  Txx 


r  5iJ00.33?f  C<f  •  Cdc^/scc2) 


(IV. E. 2-2) 


For  purposes  o£  analysis,  a  CSMP  program  (Appendix 
H)  was  written,  using  these  equations.  Figures  4.17 
through  4.24  show  the  time  responses  of  the  nz,  nz  ,  a,  B, 
q,  r,  6  ,  «v,  6„.  The  above  figures  and  table  II  show  the 

p  i  K 

time  constants  and  the  overshoots  of  the  variables. 

Figure  4.17  and  Table  II  show  that  the  achieved 
inertial  maneuver  acceleration  satisfy  the  requirements  for 
the  time  constants  and  overshoot,  when  the  cross-couplings 
which  were  discussed  before  have  been  removed. 

Figure  4.25  (taken  from  [Ref.  2]  shows  the 
critical  feedback  paths  which  couple  the  pitch  and  yaw 
channels  via  missile  roll  rate. 

The  aforementioned  figures  and  plots  verify  the 
results  of  the  [Ref.  2]. 


F.  CONTROL  LAWS  FOR  CIRCULAR  AIRFRAME 

The  control  laws  used  for  the  circular  airframe  are 

the  same  as  those  used  in  the  linear  studies,  except  for 
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‘CVI«h*  f.  *•.  •*.  r.+. oat'.' 


Fig.  4.19  Sideslip  angle  (s)  vs  Tim< 
CBTT  of  Elliptical  airfrai 
airframe  dynamics;  2  Gees 


Pig.  4.21  Body  Yaw  angular  rate  (r)  vs  T 
CBTT  Elliptical  airframe;  Idea 
airframe  dynamics;  2  Gees  (0° 
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T(SEC) 


incidence  (6  )  vs  Time  (t) 
liptical  airframe  Ideal 
ynamics;  2  Gees  (0°,  180°) 


Fig.  4.24  Roll  tail  incidence  (Gy)  vs  Time  (t) 

CBTT  of  Elliptical  airframe ;QIdealo 
airframe  dynamics;  2  Gees  (0  t  180  ) 


•.  v-v^ J. 


TABLE  IX.  CBTT  PERFORMANCE  OF  ELLIPTICAL  AIRFRAME 


VARIABLE 

VALUE  TIME  (sec) 

CROSS - 
COUPLING 

oJ 

(sec) 

ershoott 

X,  ( sec) 
Overshoott 

nZmin 

-2.934(gees) 

5.061 

All 

0.48  and 

0.41  and 

nZmax 

1.351 (gees) 

5.661 

All 

0t 

11.25% 

nZmln 

-2. 934 (gees) 

5.062 

Ideal 

0.48  and 

0.37  and 

nZmax 

0.937 (gees) 

5.853 

ideal 

0% 

0% 

amln 

1 . 35 (deg) 

6.10 

All 

amax 

6. 75 (deg) 

3.10 

All 

amln 

2. 22 (deg) 

6.295 

Ideal 

amax 

6. 75 (deg) 

3.10 

Ideal 

Bmin 

-0. 9 ( deg) 

5.49 

All 

Bmax 

1.49  (deg) 

6.05 

All 

Bmln 

1. 53 (deg) 

5.45 

ideal 

Bmax 

1.07  (deg) 

5.91 

Ideal 

^min 

-10.6(deg/sec) 

5.639 

All 

^max 

9. 5 (deg /sec) 

2.337 

All 

^mln 

-7.8(deg/sec) 

5.328 

ideal 

Snax 

9 . 5 (deg/sec ) 

2.337 

Ideal 

Pmln 

-14. 115 (deg/sec) 

6.600 

All 

Pmax 

282.094 (deg/sec) 

5.387 

All 

Pmln 

-14.291  deg/sec) 

6.600 

Ideal 

Pmax 

281.928 (deg/sec) 

5.387 

ideal 

rmln 

-2. 32(deg/sec) 

5.90 

All 

rmax 

33.1 (deg/sec) 

5. 32 

All 

rmln 

-0. 178 (deg/sec) 

7.036 

ideal 

rmax 

37. 5 (deg/sec) 

5.29 

ideal 

Note  s 

All  the 

above  referred  values  are 

for  command 

2(0°,  180°) 

1  s  62  percent  tine  constant  of  achieved  maneuver  plane 
1  acceleration  due  to  the  first  guidance  command. 

Y,:  63  percent  time  constant  of  achieved  maneuver  plane 
12  acceleration  due  to  the  second  guidance  command. 


the  gain  (a)  of  the  coordination  branch  of  lateral  control 
laws,  as  mentioned  in  the  study  for  the  elliptical 
airframe . 

The  same  2  gees  (0°,  180°)  guidance  commands  and 
flight  condition  were  used  for  direct  comparison  with  the 
performance  results  of  the  elliptical  airframe. 

However,  to  decrease  simulation  run  time,  the  first 
guidance  command  and  anti-gravity  bias  were  applied  at  zero 
time  with  no  missile  or  autopilot  initial  conditions. 

1.  Pitch  Control  Law 

The  pitch  control  law  for  the  circular  airframe  is 
the  same  as  that  use  in  the  uncoupled  autopilot  shown  in 
Figure  4.26.  The  command  used  in  the  autopilot  is 
different  than  the  one  used  in  uncoupled  channel  and  it  is 
as  follows: 


%c-  cose 

as  referred  to  in  section  IV.B.l 


(IV. B. 1-1) 


The  equations  which  represent  the  nonlinear  pitch 

control  law  are  the  same  as  those  in  section  III.C.  of 

the  linear  pitch  control  law,  using  K.  =  -0.0387 

A2 

The  equations  are: 

X  =-150  X  +•  150  (IV. F.  1-1 ) 

V  %  -1.5531  •  (v|c-0()  (iv. f. 1-2) 

Yq  *  -5^4  +  6.5?I?*(7)c-X)  (IV.  F .  1-3 ) 

<Jptz-O.W3.ypt-O.31058(y3+q/5?.3)-2.23O8(VJ-ql»3) 

(IV. F. 1-4) 
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Sf  =  -i.28.nSp  +  10*95.3  .4  (IV. f .  1-5) 

The  state  variables  of  the  above  equation  are 
shown  in  Figure  4.26. 

2.  Lateral  Control  Laws 

The  lateral  control  laws  for  the  circular  airframe 
are  the  same  as  those  used  for  the  linear  studies,  except 
the  modification  in  ag  (or  a)  which  was  referred  to  before. 

The  same  command  is  used,  as  in  the  elliptical 
airframe.  The  equations  which  represent  the  lateral 
control  laws  are  the  same  those  of  section  III.C.l  and 
III.C.2.  and  they  are: 

V2  =  -5Y2  +  1.6  Hr8  (iv. f . 2—1 ) 

Jye“  aW  -(OfP+OeP)] 


-  0.4*5  [y2+ -0.000133 .000*39  *«.p) 


(IV. F. 2-2) 


where : 


0k-  OC  7  e 

.  >  *>* 
at  -  a  ) 

a^r  ?  a<i° 

at-oS 

Sf  z  -M.H  Si  +  10*95.3  <5VC  (iv. f. 2-3) 

X,  z  -8X1  -1W  <►  +  !?.*  4>C  (IV. F. 2-4) 

Z  -5  ^  +  0.05033  (X,-P/sT.3)-hO.?5 5  C X| -P/5T -3) 

(IV. F. 2-5) 


and 


(IV. F. 2-6) 


x2  r  -6X2  +  0-0*8  (  P/5?.V) 

-15^+0.25  (VX2)+<5(V<-X2')  (iv.f.2-7) 

=  -1W^(2C  +  40*95.3  •£&.  (IV. F. 2-8) 

The  above  state  variables  are  shown  in  Figure 

4.27. 

G.  AERODYNAMIC  MODELS  FOR  CIRCULAR  AIRFRAME 

The  aerodynamic  models  of  Pitch  and  Roll-Yaw  channels 
for  the  circular  airframe  are  the  same  as  those  of  the 
elliptical  airframe  and  are  shown  in  Figures  4.1,  4.2,  the 
only  differences  which  exist  between  them,  are  the 
aerodynamic  mass  parameters  and  nonlinear  derivatives 
changes,  as  referred  to  in  the  (Appendices  A  and  B) . 

Hence,  the  equations  which  represent  the  dynamic 
models  are: 

1.  Pitch  Channel 

'Vlig  :r  -  q C* !  w  (gWs> 

q  -_|T_  +  - (dcg/sec2) 

51*3  ill 

<i  =  q  +  vc.njg-  (dtglsec) 

Substituting  the  values  for  aerodynamic  and  mass 
parameter  for  circular  airframe,  the  above  equations 
become : 

(IV. G. 1-1) 


(IV. C. 1-1) 
(IV. C. 1-2) 

(IV.C.1-3) 


-  -2.05292  •  Ctf 
9  =-0.0W52  p.r+’?36.860^CiH 


(IV. G. 1-2 ) 


Actuator 


27  CBTT  nonlinear  lateral  control  law 


(IV. G. 1-3) 


a  -  q  -  0.9$ 5H  CN-  0.01*952  p*£ 

2.  Lateral  Channels 

CCy<.€  +  Cyjy  *  <5y")  /  W  ($««s)  (IV.C.2-1) 

«  =  K.VJye  “r  (JeglseO  (IV. C. 2-2) 

V  =  ~  -gj,;-  9-  ■  q (  Cii£- €  +  Cn£y.$Y  + 

M  (IV. C. 2-3) 

P  -  V,±$*±  „  (  c<,  •  €  4  Cy^.  Sy*  C€f(t  .Sft} 

(IV. C. 2-4) 

c^d  z:  p  (deg($ec^  (iv.c.2-5) 


Substituting  the  values  in  the  above  equations, 
they  yield: 

*N&  “  2.05292  (  <*Sy  -fy)  (iv.o.2-1) 

S  z  o.qssv(ck.4+GSr-£y)  +0.01W5  pa-r  (iv.G.2-2) 
Y  ■  r -O.Oi*452  ^  P  +  *33.3*5*  (Cn*-*  +  C%*S*Y  +^.^) 

(IV. G. 2-3) 


H.  ANALYSIS  OF  NONLINEAR  CBTT  AUTOPILOTS  FOR  CIRCULAR 

AIRFRAME 

For  purposes  of  analysis,  a  CSMP  program  was  written 
(Appendix  H)  with  the  equations  representing  the 
aerodynamic  models  and  control  laws  of  the  nonlinear  CBTT 
autopilots  for  referred  airframe. 

At  first,  in  order  to  understand  the  effect  of  faster 
response,  in  the  reduction  of  the  effect  of  kinematic  and 
inertial  cross-coupling,  the  pitch  acceleration  error  gain 


0.0274  for  slower 


Ka  Figure  4.26  was  reduced  to  = 
response. 

The  equation,  which  represents  the  acceleration 


compensator  after  the  above  reduction,  is: 

V3- - — - : - hc-x) 


s-  C-S-+0 


(IV.H-l) 


Utilizing  method  of  state-space  representation  with 
forcing  term  derivatives,  the  equation  IV.H-l  yields: 

V3  =  Vit-0.9S80M-(itc-X)  (iv.H-2) 

%z-5.V  +t\.i53(\-X)  (iv.h-3) 

Using  the  above  equations,  instead  of  the  equation 
IV. F. 1-2  and  IV. F. 1-3,  a  new  CSMP  program  (Appendix  I)  was 
written. 

Figure  4.28  shows  that  the  achieved  maneuver  plane 
acceleration  has  a  good  response  to  the  first  command  which 
is  applied  at  zero  time.  The  missile  with  the  circular 
airframe  moves  upward  like  a  skid-to-turn  missile,  as  did 
the  elliptical  airframe  missile,  because  the  motion  is  in 
the  desired  maneuver  direction  and  therefore  the  roll 
channel  is  not  commanded • 

The  time  constant  of  the  achieved  maneuver  response 
is  equal  to  0.49  seconds.  The  achieved  maneuver  plane 
response  due  to  the  second  guidance  command  applied  at  3 
seconds,  shown  in  Figure  4.28  is  reacting  differently  to 
the  kinematic  and  inertial  cross-coupling  than  the 
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.  4.28  Achieved  Inertial  Acceleration  (n2)  vs  Time  (t) 

CBTT  of  Circular  airframe;  2  Gees  (0°,  180°) , 

K.  =  -0.0274 
Ai 
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elliptical  airframe  in  Figure  4.5.  Rather  than  overshoots 

and  undershoots,  a  slowing  transient  starts  at  3.5  seconds. 

Figure  4.29  shows  that  the  overshoot  in  the  body-fixed 

pitch  acceleration,  due  to  the  kinematic  and  inertial 

cross-coupling  during  the  second  command,  is  60.8  percent 

which  is  substantially  more  than  it  was  for  the  elliptical 

airframe  and  occurs  much  sooner  at  3.66  second. 

Figure  4.29  through  4.39  show  the  time  responses  of 

the  nY  ,  a,  8,  q,  p,  r,  4,  «p,  «  and  Above  figures 

c 

and  table  III  show  the  characteristic  of  the  responses  of 
the  variables  for  the  circular  ?irframe,  when  the  = 
-0.0274. 

The  referred  results  of  this  study  verify  the  results 
of  [Ref.  2] . 

I.  EFFECT  OF  INCREASING  PITCH  CHANNEL  SPEED  OF  RESPONSE 
In  order  to  reduce  the  effect  of  kinematic  and 

inertial  cross-coupling  during  the  second  guidance  command, 
the  response  of  the  Pitch  channel  of  the  CBTT  autopilot  for 
the  circular  airframe  was  made  faster  as  shown  in  section 

II.  This  was  accomplished  by  increasing 

the  acceleration  error  gain  K.  =  -0.0387. 

"2 

Using  the  above  gain  K.  the  equations  represent  the 

2 

acceleration  compensator  are  the  IV. F. 1-2  and  IV. F. 1-3. 

The  effect  of  the  change  in  KA  on  achieved  body  fixed 
acceleration  (Figure  4.40)  results  in  the  achieved  maneuver 
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Fig.  4.32  Sideslip  angle  (0)  vs  Time  (t) 
CBTT  of  Circular  airframe; 

2  Gees  (0°,  180°),  K.  =  -0.0274 
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Fig.  4.33  Pitch  Body  angular  rate  (q)  vs  Time 
CBTT  of  Circular  airframe;  2  Gees 
(0°,  180°) ,  K*  »  -0.0274 
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Fig.  4.34  Yaw  Body  angular  rate  (r)  vs  Time  (t) 
CBTT  of  Circular  airframe; 

2  Gees  (0°,  180°),  K.  =  -0.0274 
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Fig.  4.36  Roll  angle  ($)  vs  Time  (t) 
CBTT  of  Circular  airframe 
2  Gees  (0°,  180°) .  K»  *  - 


IS) 


<U>  0.1  1.0  l.S 


THCJSEC) 


Fig. 


4.38  Yaw  tail  incidence  (fiY)  vs  Time  (t) 
CBTT  Circular  airframe; 

2  Gees  (0°,  180°),  K«  =-0.0274 
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Fig.  4.39  Roll  tail  incidence  (6R)  vs 
CBTT  Circular  airframe; 

2  Gees  (0°,  180°) ,  K.  — 
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TABLE  til.  CBTT  PERFORMANCE  OF  CIRCULAR  AIRFRAME 


VARIABLE 

VALUE 

TIME (aec) 

CONDITION 

Tt  ’ sec ) I 
overshoot* 

nZaax 

1.092(gces) 

5.10 

ka 

1 

0.46  and 

nZnin 

-3. 044 (gets) 

3.060 

* 

nzaax 

1.109 (gees) 

4.540 

ka 

2 

0.34  and 

nZmln 

-3. 125 (gees) 

0.825 

5* 

nZmax 

1.002 (gaas) 

5.040 

ka 

3 

0.34  and 

nZmln 

-3. 125 (gees) 

0.825 

5* 

amax 

11.15 (deg ) 

1.432 

*A 

amin 

0 (deg) 

0 

a 

nax 

11 . 489 (deg) 

0.862 

ka 

2 

amln 

0.053(deg) 

3.656 

amax 

11. 488 (deg) 

0.862 

ka 

3 

a«iln 

0(deg) 

0 

B 

max 

4.25 (deg) 

3.461 

ka 

Bmln 

-0.399 (deg) 

4.286 

Bmax 

3.633 (deg) 

3.458 

ka 

2 

Bmin 

-0 . 326 (deg) 

4.290 

Bmax 

5.114 (deg) 

3.608 

ka 

3 

Bmln 

-0. 414 ( deg ) 

4.620 

qmax 

21. 39 (deg/sec) 

0.285 

ka 

'’min 

-8.84(deg/sec) 

3.220 

'’max 

29.562 (deg/sec) 

0.261 

"a 

2 

^ln 

-14. 323(deg/sec) 

3.225 

'’max 

29.562  (deg/aec) 

0.261 

PA 

3 

‘’min 

-17. 575 (deg/sec) 

3.251 

Paax 

278.5(deg/aec) 

3.394 

ka 

1 

Paln 

-14. 166 (deg/aec) 

4.6 

Paax 

279.487 (deg/aec) 

3.4 

«A 

2 

Pmln 

-14.071 (deg/aec) 

4.6 

P-.x 

277.574 (deg/aec) 

3.4 

Pm  in 

-13.798 (deg/aec) 

4.597 

Taax 

22.66(deg/sec) 

3.28 

rmln 

- 1 . 456 ( deg/aec ) 

4.56 

rmax 

21.486 (deg/aec) 

3.262 

rmln 

-1. 322(deg/aec) 

4.537 

rmax 

23. 159 (deg/aec) 

3.287 

rmin 

-1.364(deg/aec) 

4.818 

(sec)  ; 
overshoot* 

0.52  and 

t 

0.62  and 
0 . 05% 

0.53  and 

0t 
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plane  acceleration  response  (Figure  4.41).  Figure  4.41 
shows  that  the  achieved  maneuver  plane  acceleration 
response  during  the  first  guidance  command  has  improved. 

Table  III  shows  that  increasing  the  KA  resulted  in  a 
decreasing  in  maximum  sideslip  angle  3.66  and  a  slight 
decrease  in  yaw  angular  rates  (maximum  21.4  degrees) . 


J.  EFFECT  OF  INERTIAL  AND  KINEMATIC  CROSS-COUPLING  IN 
PITCH  CHANNEL 

The  cross-couplings,  ( —  SP )  into  a  and  (rP)  into  q  , 
were  removed  in  order  to  assess  their  effect  on 
performance.  Hence,  the  equations  for  a  and  q  in  the  pitch 
aerodynamic  model  becomes: 

a  r  q  -+  K.  Iflig  (<leg  I  sec)  (iv.j-n 

9  -5?*3q£d(M  llyy  (IV.J-2) 

Replacing  the  values  to  the  parameters  of  the  above 
equations,  they  yield: 

a  =  q  -0*9254  C #  (iv.j-3) 

9  —  ?32.260  (iv.j-4) 

where  C.,»  C  nonlinear  functions  which  vary  with  a  and 
N  m 

4  referred  to  in  (Appendix  B) .  The  lateral  aerodynamic 

r 

cross-coupling  was  retained.  For  purpose  of  analysis,  a 
CSMP  program  (Appendix  I)  was  written. 

Comparing  Figures  4.42  and  4.43  with  figures  4.40  and 
4.41  it  is  concluded  that  the  undersirable  transients  are 


gone 
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Fig.  4.41  Achieved  body-fixed  Acceleration  (n7  )  vs  Time 

(t)  CBTT  of  Circular  airframe;  Faster  respondi 
Pitch  channel;  2  Gees  (0°,  180°)  K.  =  -0.0387 
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Fig.  4.43  Achieved  Body-fixed  Acceleration  (n_  )  vs  Time 

ZB 

(t);  CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0°,  180°) , 
K.  *  -0.0387 
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Table  III  shows  that  the  achieved  maneuver  plane 
acceleration  time  constant  for  the  second  guidance  command 
has  decreased. 

Figures  4.44  through  4.5  0  show  the  time  responses  for 
a,  B,  q  ,r,  Sp*  $y'  an<^  6R  when  the  ka  *s  e<Iuai  ka=  “0*0387 
and  the  kinematic  and  inertial  cross-coupling  of  the  Pitch 
channel  have  removed. 

K.  CONCLUSIONS 

1.  The  maneuver  plane  responses  have  transients, 
which  may  have  to  be  reduced,  are  caused  by  inertial  and 
kinematic  coupling  between  pitch  and  yaw  dynamics. 

nr1  Transients  cause  excessive  overshoots  and  undershoots  in 
achieved  maneuver  plane  acceleration  for  elliptical 
airframe  (Figure  4.6)  and  excessive  slowdown  in  the  speed 
of  response  for  circular  airframe  (Figure  4.28,  4.40). 

2.  The  result  of  the  nonlinear  3-D  performance  study 
verify  the  linear  study  of  chapter  III. 

3.  The  results  of  this  study  verify  the  results  of 
[Ref.  2]. 

Further  analysis  and  conclusion  about  a,  q,  r,  p 
and  e  are  discussed  in  [Ref.  2]. 
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4.44  Achieved  Body-fixed  Acceleration  (nY  )  vs  Time 

(t)  CBTT  of  Circular  airframe;  Faster  pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  pitch  channel  removed;  2  Gees  (0°,  180°) , 
K.  =  -0.0387 
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Fig.  4.45 


Angle  of  Attack  (a)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-co 
into  Pitch  channel  removed;  2  Gees  (0°, 
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Fig.  4.46  Sideslip  angle  (0)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross¬ 
coupling  into  Pitch  channel  removed; 

2  Gees  (0,180°),  K.  =  -0.0387 
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Fig.  4.47  Pitch  Body  angular  rate  (q)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0°,  180°) , 
K*  =  -0.0387 
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Fig.  4.48  Yaw  Body  angular  rate  (r)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0°,  180°)  , 
K.  =  -0.0387 
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Fig.  4.49  Pitch  tail  incidence  ( 5p)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel) inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0°,  180°) 
K*  =  -0.0387 
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(  Sy)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0  ,  180  ) , 
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V.  MINIMIZATION  OF  THE  KINEMATIC  AND  INERTIAL 
CROSS-COUPLING  EFFECTS 


A.  GENERAL 

The  nonlinear  CBTT  autopilot  for  both  airframes, 
consisted  of  the  uncoupled  autopilots  discussed  in  section 
II,  which  coupled  via  inertial,  kinematic  and  aerodynamic 
couplings . 

The  effects  of  the  cross-coupling  were  discussed  for 
both  airframes  in  section  IV. 

In  this  section  the  kinematic,  inertial  and 
aerodynamic  coupling  effects  will  be  sumarized.  The 
kinematic  and  inertial  effects  in  the  pitch  channel  will  be 
analyzed . 

The  minimization  of  the  above  coupling  and  their 
effects  will  be  studied,  using  feedbacks  of  angle  of  attack 
and  rate  of  angle  of  attacks. 


B.  CROSS-COUPLING  OF  CBTT  AUTOPILOTS 

As  discussed  in  previous  sections,  the  four 
cross-couplings,  between  the  autopilot  channels  for  both 
airframes  are: 

1 .  Aerodynamic  Cross-Coupling 

The  aerodynamic  cross-coupling  is  the  coupling 
between  the  Roll  and  Yaw  nonlinear  channels,  via  the 
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The  major  effect  of  the 


C,  ,  C.  and  Cn  (Figure  4.2). 


above  coupling  in  the  elliptical  airframe  response  is  to 

decrease  the  overshoot  in  the  achieved  maneuver  plane 

acceleration  resulting  from  the  second  guidance  command, 

(i.e  2  gee  (180°),  the  other  effects  on  the  missile 

variables  were  small  but  in  a  direction  which  improves 

performance  (i.e  decreased  missile  body  angular  rates,  less 

control  surface  motion,  less  sideslip  variations) .  in  the 

circular  airframe,  the  effects  of  the  C,  ,  C  are 

0  6r 

negligible,  the  critical  coupling  coefficient  is  the  C„ 

6Y 

which  effects  the  stability  in  the  coordination  branch, 
which  can  be  improved  by  reducing  the  effect  of  c.  on 

S 

missile  roll  angular  acceleration.  (equation  IV. C. 2-4) 

2.  Cross-Coupling  Between  Roll-Yaw  Autopilots  for 
Coordination  Motion 

This  cross-coupling  branch,  shown  in  Figure  4.4, 
has  been  added  between  the  yaw  and  roll  control  laws,  in 
order  to  provide  coordinated  motion  between  these  two 
channels.  For  purposes  of  stability  the  coordinating 
branch  gain  Kyp  was  selected  for  the  elliptical  airframe 
Kyp  *  1  and  for  the  circular  airframe  Kyp  =  0.458.  The 
elliptical  airframe  is  more  stable  than  the  circular  with 


the  above  selected  Kyp. 
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3.  Kinematic  and  Inertial  Cross-Coupling 

The  kinematic  and  inertial  cross-coupling  exist 
between  the  pitch  and  lateral  nonlinear  aerodynamic  models 
of  the  CBTT  autopilots.  There  exist  the  following 
cross-coupling  paths  between  pitch  and  lateral  channel. 

a.  Kinematic  cross-coupling  path  of  (-Bp)  into  a 
in  the  pitch  aerodynamic  model  (Figure  4.3). 

b.  Kinematic  cross-coupling  path  of  (ap)  into  B 
in  the  lateral  aerodynamic  model  (Figure  4.4). 

c.  Inertial  cross-coupling  path  of  (pr)  into  q  in 
the  lateral  dynamic  model  (Figure  4.4). 

d.  inertial  cross-coupling  path  of(-qp)  into  r  in 
the  lateral  dynamic  model  (Figure  4.4). 

The  efffects  of  the  inertial  and  kinematic 
coupling  are  the  transients  in  the  achieved  maneuver  plane 
acceleration  response,  and  have  the  form  of  overshoots  and 
undershoots  discussed  in  section  III,  IV.  Transients  may 
also  result  in  a  slower  maneuver  plane  acceleration 
response  and  are  less  pronounced  at  higher  acceleration 
levels. 

Coupling  transients  may  be  reduced  by  increasing 
pitch  stability,  by  decreasing  maximum  roll  rate  (section 
IV),  by  improving  the  autopilot  coordination  technique  to 
minimize  sideslip  rate  8  (section  IV) ,  by  increasing  the 
effects  of  stabilizing  lateral  aerodynamic  coupling 


189 


(negative  C.  ) ,  by  increasing  pitch  channel  speed  of 
8 

response  (section  IV)  and  by  increasing  the  pitch  stability 

C  with  feedbacks  of  both  angle-of-attack  and  rate 
ma 

of  angle-of-attack  for  the  pitch  autopilot. 


C.  KINEMATIC  AND  INERTIAL  CROSS-COUPLING  IN  PITCH 
AUTOPILOT 


The  kinematic  and  inertial  cross-couplings  in  the 
pitch  autopilot,  are  via  the  pitch  angular  rate  q  and  the 
rate  of  angle-of-attack  (a)  (Figure  4.1)). 

The  equations  for  q  and  a  (figure  27)  are: 


q— Ji±5iAA 
Ivy 


r.p 

5*3 


(V.C.l) 


a  :  q  -  (P.«/5T-3) 


(V.C.2) 


where : 

rP/57.3  =  Inertia  Cross-Coupling 

-BP/57.3  *  Kinematic  Cross-Coupling 

Removing  Cm  from  the  parenthesis  in  (V.C.l)  and 
a 

combining  with  (V.C.2)  to  eliminate  q,  the  (V.C.l)  results 


in : 


a 


_  5*3 
lyv 

.  5T.35&K«»a 
Iyv 


(V.C.3) 
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Taking  the  Laplace  transform  of  (V.C.3)  and  solving 
for  a,  the  equation  becomes: 

P.8 


_  -1  .  f 

if+i  l  k,  / 


(V.C.4) 


K* 


Increasing  the  magnitude  of  Cm  will  decrease  both  K~ 

d 

and  1/K^  resulting  in  the  reduction  of  cross-coupling 
effects  on  the  angle-of-attack  (a)  and  the  influence  of  $ 
in  the  a. 


In  the  section  III  the  equation  of  n„  was  derived: 

ZB 


(V.C.5) 


(-q^/w)  ‘  ((n^*£p4  Crfa-a) 
which  involve  (a).  Therefore,  the  reduction  of  inertial 
and  kinematic  effects  on  (a) ,  results  in  the  reduction  of 


the  effects  in  the  achieved  body-fixed  acceleration  (n7  ) . 

ZB 


In  the  nonlinear  studies  the  equation  for  the  achieved 
inertial  acceleration  was  derived  as: 


%  ~  •Ca3<fr  +  ‘WyB-Sl,n4>  (V.C.6) 

It  is  apparent  from  the  equations  (V.C.4) ,  (V.C.5)  and 


(V.C.6) ,  that  an  increasing  of  Cm  results  in  the  reduction 

a 


of  kinematic  and  inertial  effects  (i.e  transients  in  the 


response)  in  the  achieved  inertial  acceleration  (n7) ,  which 


> 


it  is  desired  in  order  to  meet  the  requirement  of  time 
constants  (T^,  t2)  an^  percents  of  overshoots  and 
undershoot,  as  referred  in  Appendix  C. 

D.  MINIMIZATION  OF  THE  KINEMATIC  AND  INERTIAL 

CROSS-COUPLING  EFFECTS  IN  ELLIPTICAL  AIRFRAME 

As  discussed  in  the  nonlinear  studies,  using  the  pitch 
angular  rate  q  as  feedback  and  under  the  influence  of  the 
inertial  and  kinematic  effects,  the  achieved  inertial 
acceleration  (nz)  for  both  airframes  doesn't  meet  the 
requirements  of  the  time  constants,  overshoot  and 
undershoot,  denoted  for  the  linear  autopilots,  referred  to 
in  Appendix  C. 

Figure  4.5  and  Table  II  show  that  the  achieved 
inertial  acceleration  (nz)  for  elliptical  airframe  with 
feedback  the  pitch  angular  rate  (q)  has  overshoot  11.25% 
and  x 2  =  0.48  sec,  =  0.41  sec. 

Figures  4.40  and  Table  III  illustrate  that  the 
achieved  inertial  acceleration  ( nz )  for  circular  airframe, 
with  feedback  the  (q)  has  overshoot  5%  and  t=-.46  sec, 
t2  =  -.62  sec.  Hence  tne  (nz)  for  both  airframes  is  not 
acceptable . 

In  order  to  minimize  the  inertial  and  kinematic 
cross-coupling  effect,  namely  to  have  smoother  transient 
response  of  nz  which  will  meet  the  requirements  of  time 
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constants  and  overshoot,  the  magnitude  of  C  was  increased, 

ma 

utilizing  as  feedback  the  angle-of-attack  and  its 
derivative. 

In  order  to  increase  the  magnitude  of  C_  the 

a 

stability  augmentation  theory  was  used. 

Considered  as  stability  derivative  augmentation  o— ►Sp 
and  a—»Sp  increase  the  magnitude  of  Ma  and  M£  directly,  as 
referred  in  [Ref.  6],  the  increases  in  the  Ma  is 

9 

essentially  an  increase  in  Cm  and  hence  the  inertial  and 

a 

kinematic  effects  on  the  angle-of-attack  are  reduced. 

In  the  pitch  control  law  (Figure  3.1)  a  SAS  exists 
using  as  feedback  the  pitch  angular  rate  (q) .  In  order  to 
design  the  SAS  with  the  angle  of  attack  (a)  as  feedback, 
the  root  locus  and  lead-lag  compensation  techniques  were 
used . 

1.  Stability  Augmentation  System  (SAS)  With  Feedback 
of  the  Pitch  Angular  Rate  (q) . 

In  Figure  5.1,  the  stability  augmentation  system 
is  shown,  with  feedback  of  pitch  angular  rate  q  where  the 
transfer  function  of  the  q/6p  for  the  uncoupled  pitch 
dynamic  model,  given  in  section  II. C.  The  closed-loop 
transfer  function  of  the  above  SAS  is: 
fc(-y+l)fo3)  H>.35  5- 

G(S)= - - - — -  (V.D.l-l) 

1>  K(i+WS*3)Ht3SS-y-?1it) 

♦  0  ( s*4  o.  J  s-  1J.2T9) 
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(V.D.l-l) 


.S  . . 


Manipulating  equation  V.D.1-1,  it  results: 

Vi  (-2tt.tt3*S*-3m3fc33  _ 

a0053  S',+li)OM.sV<J.<l<-S-li.MS+  K(-Sf«i-^.32<5-?.»«*) 


(V.D.1-2) 


The  poles  and  zeros  of  the  SAS  are  as  follows: 

a.  Px  -  -188.43 

b.  P2  =  -3.45 

c.  P3  =  3.269 

d.  P4  =  0.0 
and 


e.  Z^  *  -8.0 

f.  Z2  *  -0.186 

The  other  two  zeros  are  at  infinity. 

The  root  locus  of  the  system  is  shown  in  Figure 


(5.2)  . 

The  K  was  selecting  at  the  value  of  -3.07,  because 
at  this  value,  the  system  has  two  dominant  roots  at  Pj»2  * 
-9.212**7.606  which  dominate  the  response  of  the  system. 


having  the  following  characteristics: 

».  e  =  ta-i  t ' W =  M 


where: 


b. 

c. 

d. 


J  z  COS  9=  0.733  _ _ _ 

a>  -  UJ*  vTT-5*"  =  (rouilscc.) 

=  U.203  (rad  k*c) 

w 


c  :  damping  constant 


u  :  frequency  of  system 

u  :  natural  undamped  frequency 
n 


INRC  AXIS 


The  SAS  with  feedback  the  pitch  angular  rate  q 
(Figure  5.1)  and  K  *  -3.07,  is  stable  and  the  dominant 
roots  gives  damping  ratio  c  =  0.73  in  the  response  of  the 
system. 

To  minimize  the  coupling  effects  the  angle  of 
attack  (a)  was  selected  as  feedback  in  the  previous  SAS. 

It  is  desired  the  SAS  to  have  the  same  behavior  as  before, 
for  purposes  of  stability  and  response  of  whole  the  pitch 
channel . 

Pole  zero  cancellation  technique  is  used,  to 
obtain  the  same  poles  and  zeros,  as  in  the  case  with 
feedback  the  pitch  angular  rate  q. 

For  above  purposes,  a  compensator  was  used  in  the 
SAS,  as  shown  in  Figure  5.3. 


The  transfer  function  of  the  compensator  is: 


-0.02  s-qf.35* 


(V.D.2-1) 


The  closed  loop  transfer  function  of  the  SAS  is  as 


follows: 


/-H1.3SK-MH»Wk(4 5*3  \  /-y>*s-<|f.3ST  \ 

)  \  s  ys*+o.ws-u.2?3  / 

ufl  \/-M135?-5*?.TUg\  MililW  5*5  \  f-O02-S-^.35T  \ 

1  v^3)v-0.02-s-l|t.wj‘v  l  /\^uJVsVo.W-5-iW  ) 


(V.D.2-2) 

where  the  transfer  function  a  9  -0  Q3  S-q(.35T  is 

J  s%aSfS-u.M9 

derived  in  the  section  II. C. 3. 
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Stability  Augmentation  System  (SAS)  with  a, 
Pitch  channel;  CBTT  of  Elliptical  airframe 


Rearranging  and  manipulating  the  equation 
(V.D.2-2),  it  becomes: 

«  (5.9131  sV  1229.SH  ?+  <00233.  Wq-S  4-  <32?5.0q) _ 

~  -OOOOlS5*  02215  StiH.399  S*-<M?9 S*  S  +  K  (alOW  $J+2<«U  sVf?Jtt1  S* 3iM) 

(V.D.2-3) 

The  poles  and  zeros  of  the  system  are: 

Px  =  -1877.931 
P2  =  -207.391 
P3  *  -3.447 
P4  *  -3.271 
P5  =  0.0 

and 

Z1  *  -2067.838 
Z2  »  -8.0 
Z3  »  -0.186 

The  other  two  zeros  are  at  infinity. 

The  root  locus  of  the  system  is  shown  in  figure 

5.4. 

Selecting  the  K  ■  -3.07,  the  system  has  the  same 
dominant  roots  Px  2  =  -8.249±j7.577,  as  in  section  V.D.l. 

3.  Acceleration  Loop 

Considering  the  achieved  acceleration  loop  in  the 
Pitch  channel,  the  modification  of  the  SAS,  effects  the 
transmission  gain  of  the  above  loop. 
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In  order  to  compensate  the  effect  in  the 


transmission  path,  a  compensator  was  added  as  shown  in 
Figure  5.5. 

The  transfer  function  of  the  compensator  is  the 
inverse  of  the  V.D.2-1: 

1  _  -Q.02-S-m.35? 

Gcs)=  -2|U5?-S- ?.?m  (v.D.3-1) 

Modifying  in  this  way,  the  transmission  path  of 
the  acceleration  loop,  the  gain  of  path,  is  the  same  as  in 
the  case  of  par.  V.D.l. 

4.  Measuring  of  the  Minimization  in  the  Kinematic  and 

Inertial  Cross-Coupling  Effect 

In  order  to  measure  the  minimization  a  CSMP 
program  (Appendix  J)  was  written  using  the  equations  of  the 
nonlinear  aerodynamic  models,  the  equations  of  nonlinear 
lateral  control  law  and  the  equations  of  the  modified  pitch 
control  law. 

The  equations  of  the  modified  pitch  control  law 
are  the  following: 

a.  Acceleration  lag  network 

Xr  -150X4150^  (ii.c.4-2) 

b.  Acceleration  second  order  compensator  (Figure 

V  (-0.008).  (-&03S-m.35T) 

(utvx)  ~  (.f+0  <V"D'4'1) 


5.5) 


Utilizing  inverse  Laplace  transformation  and 


rearranging  the  above  equation,  one  has: 

S+MWIYHHM  V  =  -0.000232  (UfttrX)-a<ffM-(Amt-X)  (V.D.4-2) 
Using  space  representation  of  a  system,  in  which 
the  forcing  function  involves  derivative 


[Ref.  5,  pp  675-678],  one  obtains: 


=  YI2-0.00025?5-11c  +0.000232  X  (V.D.4-3) 

W2  z  -1.11  W  W,  “  6.186?.Wj  -0.53f2  Tlc+0^??6  X  (V.  D .  4-4 ) 
V  -  (V.D.4-5) 


(V.D.4-5) 


c.  Compensator  of  Actuator 

S?c  .  . 

(y-alsi'3)  -  5  (-0.023-m.353i  (V.D.4-5a) 

Using  the  method  for  state  representation  and 


rearranging  the  equation  V.D.4-5  yields: 

r  r  .  ..  >•>.  u  /u  m I ct 


£p  r  4K-  532.3S2X10  (Y-a|5T.3)  (v.D.4-6) 

r  -206^35 ^  -  K-  1100*  W6 X 1 0C ( v. D .  4-7 ) 
<5Tf>  z  <5pt  -  K*  2S8.*J8(Y-a|51*3)  (v.d.4-8) 

d.  Equation  of  Actuator 

«5p-r  -lS*M-£p4  10?95.3^pc  (iv.b.i-4 


(V.D.4-8) 


op~  -iM.M  dp4  10?95.3-6pc  (iv.B.i-4) 

Selecting  K*3.07,  for  the  compensator  of  SAS,  with 
feedback  (a)  the  time  constant  (i.e  time  constant  of 
achieved  maneuver  acceleration  due  to  the  first  command)  is 
0.45  sec,  the  time  constant  (for  second  command)  is  0.41 
sec  and  the  undershoot  is  20.73  percent  as  shown  in  figure 


The  above  response  of  the  inertial  achieved 
maneuver  acceleration  does  not  have,  smooth  transient. 


because  of  the  large  undershoot. 

Due  to  the  above  reason,  another  K  was  selected 
equal  to  -19.37,  for  minimization  of  the  coupling  effects. 

With  this  selected  K  the  system  has  the  following 
poles,  shown  in  figure  5.4 
Px,2  *  -98  — j 9 6 . 239 
P3  =  -1879.06 
P4  =  -8.596 
P5  =  -0.189 

With  the  same  zeros  as  section  V.D.2,  the 
dominant  poles  of  the  system  are  P1#2  and  the 
characteristics  of  the  response  are: 

a.  9  *  44.24 

b.  C  *  0.716 

c.  u  *  93.775  rad/sec 

d.  wn=  134.412  rad/sec 

It  is  shown  in  known  (Ref.  6)  that  the  u  can  be 

n 

expressed  as: 


(V.D.4-9) 
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and 


(V.D.4-10) 


where: 

y:  specific  weight  of  air  1.4 
I:  moment  of  inertia 

P:  Local  pressure 

M:  Mach  number 
A:  Area 

Comparing  the  «n  a  134.412  for  the  K  *  -19.37,  it 

is  apparent  that  it  is  much  larger  than  the  <nn  *  11.203  for 

K  *  -3.07  and  hence  according  to  the  equations  V.D.4-9  and 

V.D.4-10,  the  C  increases  essentially,  resulting  in 

a 

minimization  of  kinematic  and  inertial  cross  coupling 
effects. 

Figure  5.7  shows,  that  with  K  *  -19.37  and 
feedback  in  SAS  the  (a) ,  the  achieved  inertial  maneuver 
acceleration  has  «  -.49  sec,  ^  *  "»38  sec  and  overshoot 
about  0%. 

The  above  analysis  illustrates  that  the  selection 

of  K  »  -19.37,  instead  of  K  -  -3.07,  keeps  the  damping 

ratio  about  constant  c  *  0.71  (C  *  0.73  for  K  «-3.07)  and 

increases  the  which  essentially  increases  the  C  , 
n  “*a 

resulting  in  very  smooth  response  of  the  achieved  maneuver 
acceleration. 
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Using  the  same  particular  compensators  in  the 
Pitch  channel  shown  in  Figure  5.5  and  feedback  the  rate  of 
angle  of  attack  (a) ,  the  Pitch  channel  becomes 
uncontrolable  for  both  selection  of  K  ■  -3.07  and  K  * 
-19.37,  because  the  pitch  tail  incidence  exceeds  the  limits 
of  *10°. 

Also,  selecting  the  K  *  -19.37  in  SAS  of  Figure 
5.1  with  feedback  the  pitch  angular  rate  q  the  Pitch 
channel  becomes  uncontrolable  again,  due  to  the  exceeding 
of  the  limits  in  3p. 

Therefore,  for  minimization  of  the  kinematic  and 

inertial  effects  (i.e  improving  of  the  transient  response 

of  the  achieved  maneuver  acceleration)  for  this  particular 

airframe  and  flight  conditions,  the  Cm  was  increased, 

a 

using  the  aforementioned  compensators  in  the  pitch  control 
law  (Figure  5.5)  and  feedback  the  angle-of-attack  (a). 

Comparing  the  achieved  inertial  acceleration  (nz) 
with  feedback  (q)  in  the  Pitch  channel  (Figure  4.41)  and 
the  achieved  inertial  acceleration  (nz)  with  feedback  (a) 
and  K*  -19.37  (Figure  5.7),  it  is  concluded  that  for  these 
particular  airframes  and  flight  conditions  the  feedback  of 
angle-of-attack  (a)  is  more  desirable,  than  the  feedback  of 
(q)  and  (a) ,  because  in  the  achieved  inertial  acceleration 
has  been  reduced  the  kinematic  and  inertial  effects, 
meeting  with  this  way  the  requirement  of  the  time  constants 
and  overshoot,  denoted  in  the  linear  studies. 


E.  MINIMIZATION  OF  KINEMATIC  AND  INERTIAL  COUPLING  IN 

ELLIPTICAL  AIRFRAME 

Making  a  comparison  in  Figure  5.8,  which  shows  the 

9 

kinematic  cross-coupling  (-p8)  into  <*,  in  the  pitch  channel 
with  feedback  pitch  angular  rate  (q)  and  Figure  5.10  which 
shows  the  -p8  in  the  modified  pitch  channel  (Figure  5-5) 
with  feedback  the  angle-of-attack,  it  is  concluded  (Table 
IV)  that  the  kinematic  coupling  is  increased  from 
(-0°. 988-5°. 877)  to  (-2° . 124-7° . 38 3) ,  when  the  kinematic 
coupling  effects  are  reduced  by  the  modified  pitch  channel. 

Figure  5.9  shows  the  inertial  coupling  (pr)  into  q  in 
the  nonlinear  Pitch  channel,  with  feedback  the  pitch 
angular  rate  (q) . 

Figure  5.11  shows  the  inertial  coupling  for  the 
modified  Pitch  autopilot  with  feedback  the  angle  of  attack 
(a)  (Figure  5.5). 

Comparing  the  above  figures,  it  is  concluded  (Table 
IV)  that  the  inertial  coupling  is  increased  from 
(-1°. 692-168°. 317)  to  (-20°. 373  -  181°. 415),  when  the 
inertial  coupling  effects  are  reduced  by  the  modified  pitch 
channel  discussed  in  previous  section. 

Hence,  the  minimization  of  the  kinematic  and  inertial 
effect  does  not  mean  necessarily  minimization  of  the 
cross-couplings  themselves. 
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For  minimization  of  the  kinematic  and  inertial  cross 
coupling,  the  same  modified  nonlinear  pitch  control  law  was 
used  as  that  of  section  IV. D  (Figure  5.5),  except  the 
acceleration  compensator  which  was  removed  from  the 
transmission  path  of  the  nc  as  shown  in  Figure  5.12. 

For  purposes  of  analysis  a  CSMP  program  (Appendix  j) 
was  written,  using  the  following  equation  for  the  nonlinear 
pitch  control  law  and  the  same  equation  for  lateral  control 
law  and  dynamic  models,  as  in  section  IV. B  and  IV. C. 

X  r  -  150X  +  ISO-Hie  (IV.B.1-1) 

<1  -  ~6-V-  0.5305  *V\c  +  0.43-X  (iv.b.i-2) 

5p|  -  5p2  +  V<.  532.382X403-  (V-QC|5T.3)  (v.d.4-6) 

<%  -  “20«.«S -ip*  -  V<-‘l400.mx-|06('f-als?3)  (v.d.4-7) 

5?^  -K.  259.MS  (Y-a|5*»3)  (v.d.4-8) 

<?p  r  Jp  +  4(H95.3-5pc  (iv. b. 1-4) 

where  K  *  19.37. 

Figure  5.13  shows  the  minimized  kinematic  coupling  and 
Figure  5.14  shows  the  minimized  inertial  coupling.  Figures 
5.13,  5.14  and  table  IV  illustrate  that  the  minimization  of 
the  kinematic  and  inertial  coupling  is  264%  and  195.26% 
respectively  from  its  values  in  the  nonlinear  CBTT 
autopilot  for  elliptical  airframe.  Figure  5.15  shows  the 
achieved  inertial  acceleration  (nz) ,  with  the  Pitch  channel 
modified  for  minimization  of  kinematic  and  inertial 
coupling  (Figure  5.12).  From  this  figure  it  is  concluded 
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LEGEND 


Modified  Pitch  channel  for  minimization 
of  inertial  coupling;  CBTT  of  Elliptical 
airframe 


TABLE  IV 


KINEMATIC  AND  INERTIAL  CROSS-COUPLING  FOR 
ELLIPTICAL  AIRFRAME 


CONDITION  INERTIAL  (deg/sec2)  KINEMATIC  (deg/sec2) 


Elliptical  airframe 
with  cross¬ 

MIN 

MAX 

MIN 

MAX 

coupling 
(section  IV. C) 

-1.692 

168.847 

-0.9881 

5.877 

Elliptical  ideal 

airframe  (section 

IV. E) 

Elliptical  airframe 
with  modified  Pitch 

-1.698 

167.37 

-2.124 

7.383 

autopilot 

K  =  -19.37 
(section  V.l) 

-20.373 

181.415 

00 

00 

• 

H 

1 

8.644 

Modified  Pitch 

autopilot  K  =  -30.4 
(section  V.D.) 

-7.263 

169.371 

-1.956 

8.711 

Modified  Pitch 

autopilot  for 
minimization  of 
cross  couplings 
(section  V.E) 

K  *  -19.37 

-0.0577 

57.5322 

-0.523 

1.363 

Note:  K  is  the  gain  of  the  added  compensator  in  SAS  of  the 
Pitch  autopilot. 


that  the  kinematic  and  inertial  effects  (i.e  transient  in 
achieved  acceleration  response)  is  increased. 


F.  MINIMIZATION  OF  THE  KINEMATIC  AND  INERTIAL 
CROSS-COUPLING  EFFECTS  IN  CIRCULAR  AIRFRAME 
As  shown  in  Figures  4.28,  4.40  and  4.42  the  achieved 
maneuver  acceleration  in  the  circular  airframe  has  very 
smooth  response,  in  regards  to  the  overshoot  or  undershoot 
on  it. 

The  effect  of  the  inertial  and  kinematic 
cross-coupling  in  the  nonlinear  aerodynamic  model  is  that 
the  response  of  the  inertial  achieved  acceleration  has  a 
slowing  transient  starts  at  3.5  seconds,  as  discussed  in 
section  IV. G.  Hence,  minimization  of  the  effects  means  in 
the  circular  airframe  improvement  of  time  constant  during 
the  second  guidance  command. 

1.  Stability  Augmentation  System  (SAS)  with  Feedback 
the  q. 

Using  the  same  method  for  minimization,  as  in 
section  V.D,  the  SAS  of  the  pitch  channel  for  the  circular 
airframe  with  feedback  q  is,  as  shown  in  Figure  5.16.  The 
transfer  function  of  the  q/5p  is  given  in  section  II. C. 

The  closed-loop  transfer  function  of  the  above  SAS 


V  V  • 


(*§5+  51.3  (  -  59.10**6  -  *.<3M5) 

Kl^k>  (Awr*rt(ifa+0  (V  F  H) 

<+  R.  (z^y+l)  (sy-3)  (-S9.1QM  S -^5) _ 

^(cfe+O  (&am»  S*<l<l.33K) 

Rearranging  equation  V. P.1-1,  one  obtains: 

K  (20.3*  S3+iQf0fr3Si*+303>a»fS+J34/?.y.) _ 

~  O03J -5M+ 7.0©q  5*  +  3.W9  •  S*  +  310.39  •  S  +  ^ 33  +  lc  (-<.<«,» •  S*-  40.  M- i- 

(V.F.1-2) 

The  poles  and  zeros  of  the  SAS  are:  (1) 

Pj*  -188.438;  (2)  ?2  3  -  -0 . 07 39 t j6 . 657 ;  (3)  P4  -  -0.143; 

(4)  Zx  *  -6.55  and  (5)  Z2  *  -0.118.  The  other  two  zeros 
are  at  infinity.  The  root  locus  of  the  SAS  is  shown  in 
Figure  5.17.  The  selected  value  of  K  is  -15.6  and  the 
dominant  roots  are  Plf2  *  -l®  •  736  . 536  (Pigure  5.17). 

With  those  dominant  roots,  the  response  o*  the  system  has  « 
equals  to  35.06  degrees;  C  equals  to  0.813;  -  equals  to 

7.545  rad/sec  and  «n  equals  to  13.125  rad/ sec. 

2.  SAS  with  Feedback  of  the  a 

For  minimization  of  inertial  and  kinematic  effect 
in  the  circular  airframe,  the  same  technique  was  used  (i.e 
pole-zero  cancellation  technique)  as  in  section  V.D.2  and 
the  resulted  pitch  control  law  with  feedback  of 
angle-of-attack  (a),  is  shown  in  Figure  5.18. 

The  transfer  function  of  the  added  compensator  in 
the  SAS  is  given  by  equation  V.D.2-2. 
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(V.D.2-2) 


Circular  airframe 


1HRG  RXIS 


•  Vv 


The  closed  loop  transfer  function  of  the  referred 


before  system  is  given  by  equation  V.D.2-3: 


*2= 


/  K 

h 

-S9.1-5-6.99  \ 

* 

/ -0.039-S-59.1  \ 

^-aoMS-59.1 ) 

L  S*4.0.MfS+<Hl.33  / 

K •  (^5 5T.V  (- 5W  5 -  6.99)  (-a°39  S -  59.0 
^(■^T+l)  (-°.033-5-M.O(^+0  (s+0.H8-St4H.33) 


(V.D.2-3) 


where  the  transfer  function  a/ 6  is  given  in  the 
section  II. C. 3,  rearranging  and  minor  manipulating  the 
above  equation,  one  obtains  equation  V.D.2-4: 


k.  (2&3».£+  3<Hoa2S  S*+  20 3m  65-h  226*2. 8  )  _ 

~  -aoois-sMsUiH.^i  3 -ha Atm* s-imm-tc  (o.*t+S35.m+3s*6+Hib) 

(V.D.2-4) 

The  poles  and  zeros  of  the  modified  SAS  are, 

Px  *  -1499.623,  P2  -  -188.53,  P3,4  »  -0 . 0739  ±j6 . 656  ,  P5  * 
-0.142,  Z1  *  -1499.6,  Z2  »  -6.55,  Z3  *  -0.1182  and  all 
other  zeros  at  infinity. 

The  root  locus  of  the  system  is  shown  in  figure 
5.19,  it  is  apparent  from  above  Figure,  that  for  k*-15.6 
the  dominant  roots  are  the  same  as  those  in  section  V.F.l. 

Considering  the  transmission  of  acceleration  nc 
path,  as  in  the  elliptical  airframe,  a  compensator  was 
added  in  the  path,  for  compensation  of  the  effect  in  it, 
which  is  created  by  the  added  compensator  in  the  SAS,  as 
shown  in  Figure  5.20. 


Fig.  5.19  Root  locus  of  SAS  with  a/6 
CBTT  of  Circular  airframe 


3.  Measuring  of  the  Minimization  in  the  Kinematic  and 
Inertial  Cross-Coupling  Effect 
A  CSMP  program  (Appendix  K)  was  written  using 
equations  of  nonlinear  lateral  control  law  (section 
IV. G. 2),  of  the  aerodynamic  models  (section  IV. 4)  and  the 
following  equations  of  the  modified  nonlinear  pitch  control 
law: 

a.  Acceleration  lag  Compensator 

X  r -150X  -H50  (II. C. 4) 

b.  Acceleration  Second-Order  Compensator  (Figure 

3‘20)  /  s  \  / 

^  y  o.ftfc**’  v  v~ Q-pyK-S-SM^v  F  3_i ) 

( V  5  (-f-+1 )  •  (-  59.10S-  S  -  5) 

Utilizing  inverse  Laplace  transformation  and 
rearranging  the  above  equation,  yields 

%  +5.1183  Y +0.591  Y  =*.  0.022 (\-  X  )-».9+(vX)- 

(V.F.3-2) 

Using  state-space  representation  referred  to  in 
(Ref.  5:  pp  675-678),  one  obtains: 

Wj- Ka- 0.02288  (lie- x)  (V.F.3-3) 

W2-W3-Kv  3M.?50  (mc-X)  (V.F.3-4) 

-0.5912-  W2-5.11S->^+K^.ff3J9(vx)  (v.P.3-5) 

Yj  =  (V.F.3-6) 
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c.  Compensator  of  Actuator  m  SAS 

r  /  S  .  \  / _ _  .  . 


Jfi c  _  K.(crh)  (■59.<0gi‘S-<  9a9)  (v. F« 3-7 ) 

(V  <*/5*3)  (-0*0334  5-  59.1(H) 

Utilizing,  the  aforementioned  method  for  state 
representation  and  rearranging  the  equation  V.F.3-7,  it 


yields : 


<5p4  r  <5p2  4  K  (mS896.0JBS)*  (Y3- aisj-3)  (v.F.3-8) 
5p2  S-an**W3-^  -  1S00  ^-y.7.3X<0TfH»5)(V.  F.  3-9 ) 
5fc=  <5ft  -  K.  (32.^4)  (Y3- a/573)  (V.F.3-10) 

Equation  of  Autuator 

P.  r  r»  (IV. G. 3-8 

ip  =  -1884  *6p4*  <0795.3- ipQ 


(IV. G. 3-8) 


Selecting  the  K  *  -15.6  and  K*  *  -0.0387 

2 

(Figure  4.26)  and  feedback  the  pitch  angular  rate  (q) ,  the 
achieved  inertial  acceleration  nz ,  has  =  -.34  sec,  T2  “ 
0.62  sec  and  overshoot  5%  (Figure  4.40,  Table  III). 

This  response  of  the  nz  is  unacceptable 
because  it  doesn't  meet  the  requirement  on  the  time 
constant,  referred  to  in  Appendix  C. 

Selecting  K  *  -15.6  and  KA  =  -0.0387  (figure 
5.20)  and  feedback  the  angle  of  attack  (a)  in  the  pitch 
channel  the  achieved  inertial  acceleration  nz  (Figure  5.21) 
has  a  -.35  sec,  t  =  0.6  sec  and  overshoot  7%,  which  is 
again  unacceptable. 


(V.G.3-9) 


It  is  known  (Ref.  6)  that: 


Hence,  in  order  to  reduce  the  *2  it  is  needed 
to  increase  the  Cm  .  Selecting  K  =  -70.88  and  KA  =  -0.0387 

d 

(Figure  5.20)  with  feedback  the  angle  of  attack  (a) ,  the 
dominant  roots  of  the  SAS  (Figure  5.18)  are  p^,2  = 

-90. 55±j87.95  and  the  characteristics  of  the  system 
response  are  0  =  44.16  degrees,  ;  =  0.79  ,  «•>  =  87.932 
rad/sec  and  «>n  =  126.212  rad/sec.  With  that  selected  k  the 
z  decreases  from  0.818  (when  K«  -15.6)  to  0.72  and  the  « 
increases  from  13.125  rad/sec  (when  K  *  -15.6)  to  126.212 
rad/sec. 

Figure  5.22  shows,  the  achieved  inertial 

acceleration  n7,  when  K  *  -70.88,  K.  =  -0.0387  and  feedback 
2  A2 

the  (a).  The  f ^  »  0.33  sec,  *2  *  0.48  sec  and  overshoot 
14.33  %.  It  is  apparent  that  the  ^2  decreases,  which  is 
desirable  but  the  overshoot  increases,  because  the  damping 
ratio  c  decreases. 

In  order  to  reduce  the  overshoot  and  to  keep 

the  time  constants  below  the  0.5  sec,  the  KA  reduces  to 

-0.0274  (Figure  5.20)  and  the  gain  Kyp  (Figure  4.27)  of  the 

coordination  branch  in  the  nonlinear  lateral  control  law 

increases  from  0.458  to  1.0. 

Using  Kyp  ■  1.0  (Figure  4.27)  in  the 

coordination  branch  in  the  lateral  control  law,  K. 

A1 
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CBTT  of  Circular  airframe;  K  =  -70.88 
K.  =  -0.0387 


Achieved  Inertial 


CBTT  of  Circular  airframe;  K  =  -70.88 
K.  =  -0.0274 


*  -0.0274  in  the  acceleration  compensator  and  K  =  -70.88  in 
the  added  compensator  in  the  SAS  with  feedback  the 
angle-of-attack  (a)  (Figure  5.20),  the  obtained  achieved 
maneuver  acceleration  has  ^  =  -.43  sec,  *2  =  -.37  sec 
overshoot  8.9%,  as  shown  in  figure  5.25  and  table  V. 

Using  in  the  pitch  channel  with  the 
aforementioned  particular  compensators,  as  feedbback  the 
rate  of  angle  of  attack  (a)  and  making  all  the  above 
selection  of  ka,  K  and  Kyp,  the  whole  system  becomes 
uncontrolable  because  the  required  tail  incidence  angle 
exceeds  the  limits  of  ±10°. 

Also  using  the  selections  of  k,  Ka,  Kyp  and 
feedback  pitch  angular  rate  q  (Figures  4.26,  4.27),  again 
the  whole  system  becomes  uncontrolable,  because  the 
required  tail  incidence  angle  exceeds  the  limit  of  ±10°. 

The  above  analysis  illustrates  that  the 
angle-of-attack  is  more  desirable  feedback  than  the  (q)  and 
(a) ,  in  the  nonlinear  pitch  control  law  for  these 
particular  circular  airframe  and  flight  conditions,  because 
it  reduces  the  kinematic  and  inertial  coupling  effect,  such 
that  the  achieved  inertial  acceleration  (nz)  meets  the 
requirements  of  the  time  constants  and  overshoot,  referred 
to  in  the  linear  studies. 
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TABLE  V 


CHARACTERISTIC  OF  THE  INERTIAL  ACHIEVED  ACCELERATION  (nz) 
FOR  CIRCULAR  AIRFRAME 


CONDITION 


(sec)  tj  (sec)  overshoot 


ka  31 

0.0387, 

K  =  -15.6 

0.35 

0.6 

7 

ka  38 

0.0387, 

K  *  -70.88 

0.33 

0.  48 

14.33 

ka  * 

0.0274, 

K  =  -15.6 

0.435 

0.51 

8.6 

ka  = 

0.0274, 

K  =  -70.88 

«vd= 

1 

0.43 

0.37 

8.9 

Note : 


K.:  pitch  autopilot  acceleration  gain  in  acceleration 
A  compensator 


r  :  63  percent  time  constant  of  achieved  maneuver 
plane  acceleration  due  to  first  command 

T 


2 .  63  percent  time  constant  of  achieved  maneuver 
plane  acceleration  due  to  second  command 


KYp:  gain  of  the  coordination  branch 

K:  gain  of  the  added  compensator  is  the  SAS  of  pitch 
autopilot 
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G.  MINIMIZATION  OF  KINEMATIC  AND  INERTIAL  COUPLING  IN 
CIRCULAR  AIRFRAME 


As  in  the  elliptical  airframe  comparing  the  values  of 
the  inertial  and  kinematic  cross-coupling  for  the  circular 
airframe,  namely,  the  (P.r)  and  (-p8) ,  referred  to  in  table 
VI,  it  is  concluded  for  this  particular  analysis  that  as 
the  coupling  effects  in  the  nonlinear  Pitch  autopilot  are 
reduced,  the  cross-couplings  themselves  are  increased. 

In  order  to  minimize  the  cross-couplings,  the  pitch 
control  law  was  modified,  as  shown  in  Figure  5.26. 

For  purposes  of  analysis,  a  CSMP  program  (Appendix  K) 
was  written,  using  the  following  equations  for  the 
nonlinear  modified  pitch  control  law,  lateral  control  law 
and  dynamic  models,  as  in  sections  IV. G  and  IV. H. 


1*  i  r  H50X  +<50  <*18  (II. c. 4) 

2.  r  W2 0.022SS(>tc'^  (v.f.3-3) 

3.  W2  r  W3- Ka- 34.850 (tu-X)  (v.f.3-4) 

4  •  W3  =  -0.591  W2 -  5.  m W3  +  *A  ‘  ^3.385  (\-X)  (v .  f.  4-5 ) 

5.  V3  -  Wl  (V.F.4-6) 

6.  5pt  =  -0.ma-Spc-0.3H058(Va/s7.3)-2230^*(y3-£X/5?.3) 

(V.F.4-7) 


which  is  the  equation  IV. G. 1-4,  except  of  the 
feedback,  which  is  the  angle  of  attack  (a)  instead  of  the 
pitch  angular  rate  (q)  . 

7.  <fp  r  410*95.3- £pc  (iv.g.i-5) 
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TABLE  VI 


KINEMATIC  AND  INERTIAL  COUPLING  FOR  CIRCULAR  AIRFRAME 


CONDTION 

INERTIAL 

(deg/sec2) 

KINEMATIC 

(deg/ sec  ) 

MIN 

MAX 

MIN 

MAX 

Circular  airframe 
with  cross¬ 
couplings 
(section  IV. G) 

-0.008 

92.794 

-1.132 

0.012 

with  cross¬ 
coupling  and 
faster  response 

K  -  -0.0387, 

K*  -  -15.6 
(section  V.F) 

-0.0768 

100. 573 

-2.04 

0.0353 

with  cross¬ 
coupling  and 
faster  response 

Ka  *  -.0387, 

K  -  -70.88 
(Section  V.F) 

-0.0212 

103.553 

-1.89 

0.0436 

with  cross¬ 
coupling  and 
slower  response 

K.  *  -0.0274, 

K*  -  -15.6 
(section  V.F) 

-0.107 

118.294 

-2.97 

0.0387 

with  cross¬ 
coupling  and 
slower  response 

K  -  -0.0274, 

K*  «  -70.88 
(section  V.F) 

-0.0743 

109.411 

-2.32 

0.0431 

with  cross¬ 
coupling 
slower  response 
and  increased  gain 
in  the  coordinated 
branch  Kvp  -  1.0, 
Ka  »  -0.287 

K  »  -70.88 
(section  V.F) 

-6.0791 

199.841 

-0.8662 

1.0541 

with  cross¬ 
coupling  slower 
response  and 
modified  pitch 

-2.264 

38.386 

-0.00788 

-0.6021 

control  law  as 
saction  V.G 
Ka  «  -0.0274 


Note:  Ka,  K,  Ky„  as  referred  in  section  V.F  and  V.G. 
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Figures  5.27  and  5.28  show  the  minimized  inertial  and 
kinematic  cross-coupling,  with  feedback  the 

angle-of-attack,  with  K.  *  -0.0274  and  K  =  -70.88. 

1 

Comparing  the  values  of  the  above  minimized  inertial  and 
kinematic  coupling,  with  the  values  of  the  kinematic  and 
inertial  coupling  when  the  pitch  channel  has  feedback  the 
pitch  angular  rate  (q)  and  KA  =  -0.0274,  as  given  in  table 
VI,  it  is  concluded  that  the  couplings  when  the  pitch 
channel  has  feedback  angle-of-attack  (a)  and  the 
aforementioned  particular  compensator  (Fgure  5.26),  were 
reduced  228.29%  for  inertial  and  190%  for  kinematic. 

Also,  comparing  the  Figure  5.25  which  illustrates  the 
achieved  inertial  acceleration  nz  for  circular  airframe 
with  minimized  coupling  effects  and  Figure  5.29  illustrated 
nz  for  circular  airframe  with  minimized  coupling,  it  is 
concluded  that  as  the  coupling  reduces,  the  coupling 
effects  increases. 

H.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  that  follow  are  based  on  a  single 
representative  flight  condition  (Mach  3.95  and  60  kft 
altitude)  and  for  the  particular  analysis  which  was  done  in 
this  chapter. 

1.  From  the  studies  in  this  chapter,  it  is  concluded: 
a.  Transients  in  maneuver  plane  acceleration  are 


Fig. 


5.29  Achieved  Inertial  acceleration  (n~)  vs  Time  ( t) 
Modified  pitch  channel  for  minimization  of 
coupling;  CBTT  of  Circular  airframe 
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caused  by  kinematic  and  inertial  coupling  between  pitch  and 
yaw  dynamics  through  missile  roll  rate.  Transients  are  in 
the  form  of  overshoots  and  undershoots  (section  V.B) . 

b.  Using  as  feedback  in  the  nonlinear  pitch 
channel  for  both  airframes  the  pitch  angular  rate  (q) ,  the 
transients  in  the  maneuver  plane  exceeds  the  denoted 
requirements  and  the  achieved  maneuver  acceleration  is  not 
acceptable  (sections  V.D  and  V.F) . 

c.  Using  as  feedback  in  the  nonlinear  pitch 
channel  for  both  airframe,  the  angle-of-attack  (a)  and 
modifying  suitably  the  pitch  control  law  (Figure  5.5  and 
Figure  5.20),  the  transients  are  minimized  below  the 
desired  requirement  and  the  achieved  inertial  acceleration 
(nz)  is  acceptable.  (sections  V.D  and  V.F). 

d.  Using  as  feedback  the  rate  of  angle  of  attack 
(a) ,  the  whole  system  becomes  uncontrolable  because  th« 
pitch  tail  incidence  (<$p)  exceeds  the  limits.  (sections 
V.D,  V.F) . 

e.  Using  as  feedback  the  rate  of  angle  of  attack 
(a)  and  modifying  suitably  the  pitch  control  law  (Figure 
5.26),  the  kinematic  and  inertial  coupling  are  reduced 
dramatically,  comparing  with  their  values  when  the  used 
feedback  is  the  pitch  angular  rate,  (sections  V.D,  and  V.G) 

f.  Using  as  feedback  the  rate  of  angle  of  attack 
(a),  in  the  modified  pitch  channel  (Figure  5.26)  for 


minimization  of  kinematic  and  inertial  coupling,  the  whole 
system  becomes  uncontrolable ,  because  the  pitch  tail 
incidence  ( 6^)  exceeds  the  limits,  (sections  V.E  and  V.G) 

g.  Minimizing  the  inertial  and  kinematic 
cross-coupling  effects,  the  values  of  kinematic  (-pB)  and 
inertial  (pr)  coupling  increase,  (section  V.E  and  V.G) 

h.  The  above  conclusions  show,  that  for  the 
particular  airframe,  flight  conditions  and  analysis  is  not 
possible  to  achieve  minimization  of  the  kinematic  and 
inertial  cross-coupling  effects  and  simultaneously 
minimization  of  the  cross-couplings  themselves. 

2.  Recommendations 

Considering  the  studies  of  this  work  and  the 
conclusions,  it  is  recommended: 

a.  Further  analysis  of  the  minimization  of  the 
kinematic  and  inertial  cross-coupling  effects  and  the 
cross-coupling  themselves  in  other  airframes,  and  with 
other  analysis  to  achieve  simultaneous  minimization. 

b.  Further  analysis  of  the  CBTT  autopilots  in 
other  flight  conditions  for  verification  of  the  above 
results. 
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APPENDIX  A 


MISSILE  SIZING  AND  MASS  PROPERTIES 
In  order  to  provide  a  realistic  missile  based  on  the 
configuration  concepts  tested  aerodynamically  in  [Ref.  1]  , 
the  models  have  the  following  characteristics. 


Characteristic 

Circular 

Elliptical 

Length  (in) 

168 

168 

Max.  Diameter 

24 

Max.  Major  axis  (in) 

41.57 

Max,  Minor  axis  (in) 

13.86 

Weight  (lbs) 

2525 

2475 

Xxx  (slug-ft2) 

40 

110 

Iyy  (slug-ft2) 

804 

790 

Iz2  (slug-ft2) 

810 

853 

c.g  distance  from  L.E  (in) 

100.8 

100.8 

Reference  length  d  (ft) 

2 

2 

2 

Reference  area  S(ft  ) 

TT 

IT 

APPENDIX  B 


NONLINEAR  AERODYNAMIC  DATA 


The  aerodynamic  data  were  taken  from  [Ref.  2].  For 
reference,  the  normal  and  pitching  moment  carves  at  M  = 
3.95  have  been  produced  in  fig.  B.l  and  B.2  for  the  two 
configurations  in  [Ref.  2].  The  aerodynamic  derivativies 
Cv,  C  ,  C.,  with  respect  to  sideslip  angle  g,  yaw  control 
5„ ,  and  roll  control  are  presented  in  figures  B.3,  B.4 

I  K 

and  B.5  of  [Ref.  2]  as  used  in  the  computer  simulation, 
namely,  as  piece  wise  linear  segments  for  ease  in 
interpolation. 
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APPENDIX  C 


REQUIREMENTS  FOR  UNCOUPLED  AUTOPILOTS 


The  requirements  for  the  classical  design  technique  of 
the  uncoupled  autopilots  are  given  in  [Ref.  2]  and  are  the 
following : 

1.  Acceleration  time  response 

a.  63  percent  time  constant  of  0.5  seconds  for  a 
step  command  of  acceleration  at  the  flight  condition  of 
interest,  (i.e  M  ■  3.95  and  altitude  60kft)  and  small 
angles-of-attack.  This  response  is  representative  of  a 
tactical  missile  of  this  size. 

b.  Overshoot  10  percent 

c.  Zero  steady  state  error  in  acceleration  to 
reduce  variations  of  guidance  navigation  gain. 

2.  Relative  stability 

Gain  margins  6db,  phase  margins  30  deg  with  a 
goal  of  12  db  and  50  deg. 

3.  High  Frequency  Attenuation  in  Actuator  Command 
Branch 

It  must  be  ^  15  db  at  100  rad/sec  and  zero 
angle-of-attack.  This  requirement  will  provide  sufficient 

high  frequency  attenuation  for  30  Hz  actuator  and  for  body 
bending  modes,  when  high  frequency  filters  are  added.  This 

requirement  also  limits  autopilot  speeds  of  response. 
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CSPP  PROGRAMS  FOR  PITCH  UNCOUPLED  AUTCPILOT 
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PITCH  EQUATIONS  OF  CONTROL  LAN 
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